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Preface 



Automated deduction in geometry has several roots that go back to developments 
in the 20th century. These fundamental methods in ADG are the quantifier 
elimination method of Tarski and Collins, the method of characteristic sets of 
Ritt and Wu, and the Grobner basis method of Buchberger. Based on these 
algorithmic techniques various geometric theorem provers have been developed 
in recent decades. Applications of ADG theorem provers range from computer- 
aided geometric design to robotics and education. 

In 1992 H. Hong, D.M. Wang, and F. Winkler organized the workshop “Al- 
gebraic Approaches to Geometric Reasoning” in Schloss Weinberg near Linz, 
Austria. The idea of this workshop was to bring together mathematicians, com- 
puter scientists and people working in applications, to discuss the newest results 
in ADG. After this first workshop D.M. Wang took over the responsibility of 
organizing a series of workshops in ADG. He organized and chaired ADG 1996 
in Toulouse. The following workshops were ADG 1998 in Beijing organized and 
chaired by X.S. Gao, and ADG 2000 in Zurich organized and chaired by J. 
Richter-Gebert. The proceedings of ADG 1996, ADG 1998, and ADG 2000 have 
all appeared in the series Lecture Notes in Artificial Intelligence. 

I was asked by D.M. Wang to organize and chair ADG 2002. 1 gladly accepted 
this offer, and ADG 2002 took place on September 4-6, 2002, in Schloss Hagen- 
berg, the home of the Research Institute for Symbolic Computation ( RISC-Linz ) 
of the Johannes Kepler University of Linz, Austria. The speakers at the work- 
shop were invited to revise their work and submit papers to the proceedings. 
These submissions were reviewed by the program committee and 13 papers were 
selected to be included in the proceedings of ADG 2002. Let me just briefly 
review the various contributions to these proceedings. 

Several of the papers are concerned with new theoretical developments in 
ADG: G. Bodnar investigates an important subproblem in the resolution of 
singularities, namely the question of whether two (or several) varieties have the 
normal crossing property. H. Li presents an innovative new representation of 
objects in affine and projective algebraic geometry using Cayley algebra and 
bracket algebra. J.C. Owen and S.C. Power deal with a configuration problem 
of points in a plane. Such configuration problems can be phrased in terms of 
graphs. They show that certain related graph problems cannot be solved by 
radicals. I.J. Tchoupaeva proposes a coordinate-free representation of varieties 
in Grassmann algebra along with Grobner basis techniques for computing in 
such an environment. X. Chen and D.K. Wang use a projection algorithm for 
quasi varieties in geometry theorem proving and show that their method provides 
optimal nondegeneracy conditions for geometric theorems. 

Another group of papers is concerned with software and computer systems 
for geometric theorem proving: X.S. Gao and Q. Lin describe their software pack- 
age MMP/Geometer, implementing Wu’s method for Euclidean and differential 




VI 



Preface 



geometry. H.G. Grabe presents the SymbolicData GEO records, a proposed stan- 
dard for proof schemes. A. Pasko and V. Adzhiev design a specialized language 
for shape modeling. D.M. Wang has further developed his system GEOTHER, 
and shows how it can be applied for solving geometric proving problems. 

Finally, there is a group of papers investigating applications of ADG to prob- 
lems in computer-aided geometric design and robotics: C. Jermann, B, Neveu, 
and G. Trombettoni consider the problem of structural rigidity in geometric 
constraint satisfaction problems. They propose a new approach to this concept, 
which takes into account geometric aspects of rigidity. M. Shalaby, B. Jiittler, 
and J. Schicho propose a new method for constructing a low degree G 1 implicit 
spline representation of a given parametric planar curve. L. Yang in his paper 
introduces an invariant method based on distance geometry to construct the con- 
straint equations for spatial constraint solving. S. Corvez and F. Rouillier use 
computer algebra tools to give a classification of 3-revolute-jointed manipulators 
based on the cuspidal behavior. 

I want to thank the RISC-Linz institute for its hospitality during the work- 
shop, and the Government of Upper Austria, the University of Linz, the Linzer 
Hochschulfonds, and the Wirtschaftskammer Osterreich for their financial sup- 
port. Last but not least I want to thank Manasi Athale for her work in collecting 
all these contributions and preparing them for printing. 

The series of workshops in ADG will continue, with ADG 2004 being orga- 
nized by Neil White at the University of Florida. 



October 2003 



Franz Winkler 
Johannes Kepler Universitat 
Linz, Austria 




Table of Contents 



Algorithmic Tests for the Normal Crossing Property 

Gabor Bodnar 1 

The Projection of Quasi Variety and Its Application on Geometric Theorem 
Proving and Formula Deduction 

XueFeng Chen and DingKang Wang 21 

Using Computer Algebra Tools to Classify Serial Manipulators 

Solen Corvez and Fabrice Rouillier 31 

MMP/Geometer - A Software Package for Automated Geometric Reasoning 
Xiao-Shan Gao and Qiang Lin 44 

The SymbolicData GEO Records - A Public Repository of Geometry Theorem 
Proof Schemes 

Hans-Gert Grabe 67 

A New Structural Rigidity for Geometric Constraint Systems 

Christophe Jermann, Bertrand Neveu, and Gilles Trombettoni 87 

Algebraic Representation, Elimination and Expansion in Automated Geometric 
Theorem Proving 

Hongbo Li 106 

The Nonsolvability by Radicals of Generic 3-connected Planar Graphs 
John C. Owen and Steve C. Power 124 

Function-Based Shape Modeling: Mathematical Framework and Specialized 
Language 

Alexander Pasko and Valery Adzhiev 132 

C 1 Spline Implicitization of Planar Curves 

Mohamed Shalaby, Bert Jiittler, and Josef Schicho 161 

Analysis of Geometrical Theorems in Coordinate-Free Form by Using 
Anticommutative Grobner Bases Method 

Irina Tchoupaeva 178 

GEOTHER 1.1: Handling and Proving Geometric Theorems Automatically 
Dongming Wang 194 

Distance Coordinates Used in Geometric Constraint Solving 

Lu Yang 216 



Author Index 



231 




Algorithmic Tests for the Normal Crossing 

Property 



Gabor Bodnar 

Research Institute for Symbolic Computation 
Johannes Kepler University, A-4040 Linz, Austria 

bodnar@r isc . uni-linz . ac . at 



Abstract. Given a finite set of varieties in some nonsingular affine va- 
riety W. They are normal crossing if and only if at every point of W 
there is a regular system of parameters such that each variety can be 
defined locally at the point by a subset of this parameter system. In this 
paper we present two algorithms to test this property. The first one is 
developed for hypersurfaces only, and it has a straightforward structure. 
The second copes with the general case by constructing finitely many 
regular parameter systems which are “witnesses” of the normal crossing 
of the varieties over open subsets of W. The ideas of the methods are 
applied in a computer program for resolution of singularities. 



1 Introduction 

The notion of normal crossing is ubiquitous in singularity theory and it has 
a prominent role in the problem of resolution of singularities (see e.g. [1] for 
the foundations of the classical theory, or [2] for a more recent survey on the 
field) and in constructing alterations (see e.g. [3]). Informally speaking, finitely 
many subvarieties of some nonsingular variety W are normal crossing iff at every 
point of W there exists an “algebraic coordinate system,” such that each of the 
varieties can be defined locally at the point by coordinate functions. This implies 
in particular that the varieties themselves cannot have singularities. From this 
we can obtain the intuition that the normal crossing singularities (that is, points 
where normal crossing varieties intersect) are, so to say, the best ones after 
nonsingular points. 

This is also reflected in the definition of embedded resolution of singularities, 
which requires that a desingularization morphism 7 r : W' — >• W of an embed- 
ded variety X C W (where W, W' are nonsingular) has the property that the 
irreducible components of the total transform 7r -1 (X) are normal crossing and 
the strict transform, i.e. the Zariski closure of n~ 1 (X \ S), is nonsingular. Here 
S C A is the set of singularities of X and tt is an isomorphism outside S. As 
the matter of terminology, 7 r _1 (S l ) is called the set of exceptional components, 
which is a union of hypersurfaces in W' . 

In the resolution algorithm of Villamayor (for the theoretical presentation 
see [4, 5], for the implementation see [6, 7, 8]) it is essential to keep the excep- 
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tional components normal crossing in the intermediate steps, when the desin- 
gularization map 7r is constructed as the composition of well chosen elementary 
transformations, called blowups. 

It is not difficult to choose blowups that create new exceptional components 
which are normal crossing with the old ones. On the other hand, it is not always 
possible to have such transformations that at the same time also reduce the 
multiplicities of the singularities. Therefore the algorithm often has to apply 
“preparatory” blowups which do not reduce multiplicity; whose only job is to 
create an advantageous situation, in which the next transformation can achieve 
both goals at the same time. 

The subtask of the resolution which deals with this “preparation” is usually 
referred to as the transversality problem. The simplification of the solution of this 
problem motivated the investigation of the normal crossing property. The ideas 
and parts of the algorithms of this paper are applied in a recent implementation 
of the resolution algorithm. Their integration into the theory of Villamayor is 
beyond the scope of this presentation, and is planned to be described in an 
forthcoming paper. 

In Sect. 3 we present an easy algorithm to test the normal crossing of hy- 
persurfaces, which cannot be generalized in a straightforward way to settle the 
question for varieties. Therefore in Sect. 4 we start with a different approach 
to tackle the general case. The resulting algorithm is more complicated than 
the one of Sect. 3; on the other hand, it does not only say “yes” or “no,” but 
it also computes a special representation of the input varieties. This consists of 
finitely many “algebraic coordinate systems” (i.e. regular parameter systems, see 
Sect. 2) together with corresponding open subsets and tags for each variety. The 
representation has the property that for each “coordinate system,” a given vari- 
ety is defined, over the open subset belonging to the system, by the coordinate 
functions which are marked by the tag of the variety. 

The prerequisites for the reader to follow the paper are elementary results 
from commutative algebra, algebraic geometry and differential geometry. 

To my best knowledge, no algorithmic solution has ever been given for the 
subject of this paper. I would like to thank Herwig Hauser for his helpful remarks 
and suggestions. This work has been supported by the Austrian Science Fund 
(FWF) in the frame of the project “Explicit Resolution and Related Methods 
in Algebraic Geometry and Number Theory” (P15551). 

2 Preliminary Definitions 

For the sake of simplicity we take k to be an algebraically closed field of char- 
acteristic zero (e.g. the complex numbers). Actually, the methods of the paper 
would also work for algebraically closed fields of positive characteristic for almost 
every particular example. In characteristic p the failing cases would arise from 
the collapse of partial derivatives of pth powers. 

We denote the affine m-space over k by A™, with coordinate ring R = 
k[ x \, . . . , x m }. An affine algebraic set W C A™ is defined by a polynomial ideal 
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J Q R, where, because of noetherianity, we can always take finite generating sets 
for J. The coordinate ring of W is Ow = R/J. The Krull dimension dim Gyp 
is the supremum of the height of all prime ideals in Ow, where the height of 
a prime ideal P is the supremum of all integers n for which there exists a chain 
of prime ideals To C • • • C P n = P C Ow- The dimension dim W is by definition 
dim Ow ■ 

We define the Zariski topology on an algebraic set W by taking the algebraic 
subsets of W to be the closed subsets. An algebraic set is called a variety iff it 
is irreducible in the Zariski topology. 

At a point a £ W, the local ring Ow,a is the ring of regular functions at a, 
i.e. the ring of Ow - rational functions with nonvanishing denominator at a. The 
(unique) maximal ideal of Ow.a is denoted by m a , which, by default, can be 
generated by x\ — ai, . . . , x m — a m if a = (ai, . . . , a m ). The local ring Ow.a is 
called regular if m a can be generated by dim Ow.a many elements. 

The algebraic set W is nonsingular at a point a £ IT iff Ow,a is regular, 
and W is nonsingular iff it is nonsingular at all of its points. 

Definition 1. Let Ow,a be a regular local ring, a minimal set of generators of 
m a is called a regular system of parameters of Ow, a- 

Let n = dim W, so that p = (p±, . . . ,p n ) is a regular system of parameters 
of Ow,a , where a is a nonsingular point of W. Then there exists an open neigh- 
borhood U C W of a such that for all a' £ U, pi — Pi(a '), . . . ,p n — p n (a’) is a 
regular parameter system of Ow,a' (see Remark 1). In this case we call p a regular 
parameter system over U. 

We assume from now on that W is a nonsingular variety. Then without loss 
of generality we can assume that we have regular parameter systems over W, 
otherwise we could just pass to open subsets of a suitable covering of W. 

If p is a regular system of parameters over W, we have dpi , . . . , dp n forming 
a basis of the module of differential forms fiw/k- Accordingly, they give rise to 
partial derivations d Pl ,...d Pn £ Derfc(W). To compute derivatives over W it 
is enough to know the Jacobian of the generators of Ow, that is P £ 
where Pij = d Pi Xj. Then for a representative / £ k[x i, . . . ,x m ] of an element 
from Ow, dpif = Xqli Pij^xjf represents its partial derivative. From now on, 
whenever we have a system of regular parameters over an affine variety, we 
assume to have also the matrix P available, so that we can effectively compute 
partial derivatives. 

Remark 1. Given a regular system of parameters p of Ow,a, the largest open 
subset U C W over which p is (or with more precise terminology: gives rise to) 
a parameter system consists of the points where none of the denominators in p 
vanish and the vector field defined by dp i, . . . , dp n is linearly independent. The 
latter condition follows from the fact that p is a regular system of parameters 
at a iff (the image of) p is a basis of the fc-vector space m a /m„ (by Nakayama’s 
lemma) . This is the case when the vectors defined by dpi , . . . , dp n are linearly 
independent at a. 
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Next we provide other means to compute the singularities of subvarieties 
of W, which we will use in the algorithms of the upcoming sections. Let X C W 
be a subvariety of codimension r (= n — dim X), defined by /i, . . . , fi £ Ow- 
The Jacobian matrix of X with respect to the defining equations fi, ... , fi is 

dw(/li ■ ■ • j fl) = (.^Pj fi)j=l,...,n;i=l,...,li 

and the Jacobian ideal Uiw,r(/i> • ••>/;) is generated by all the r x r minors of 
the Jacobian matrix. The singular locus Sing(X) is then defined by (/i, . . . , fi} + 
3W,r(/ij ••■>//) in Off- 

In the case of hypersurfaces the previous computation simplifies in the fol- 
lowing way: Let X C W be a hypersurface defined by / £ Ow , then Sing(X) is 
defined by (/) + (d Pj f \ j = 1, . . . , n). 

Definition 2. Let X\, ... , X s C W be varieties. They are normal crossing at 
a point a £ W iff there exists a regular parameter system p at a, such that 
each of the Xi can be defined, locally at a, by a ( possibly empty) subset of this 
parameter system. We say that Xi , . . . , X s are normal crossing ( over W ) iff they 
are normal crossing at every point (ofW). 

Remark 2. (i) If X \ , . . . , X s are normal crossing at a point a, then they are 
normal crossing in an open neighborhood of a. In other words, the “is normal 
crossing” predicate, considered as a function, is upper semi-continuous, that is, it 
can jump from true to false only at special points. This property follows from the 
remark after Definition 1. (ii) Being nonsingular is a necessary condition for all 
the varieties in order to be normal crossing, because a variety at a singular point 
cannot be defined by regular parameters, (iii) The definition can be extended 
to algebraic sets or even to subschemes of W. First of all, we would pass to 
the corresponding reduced subschemes, then to their irreducible components to 
reduce the problem to the one we gave above. This is the reason why we prefer 
to use the simpler language of varieties in this paper. 

In Fig. 1, we show a few negative and positive examples of the normal cross- 
ing property. The points of “non-normal crossing” are marked with black. In 
Fig. 1(a) all the singular points are not normal crossing. The surfaces in Fig. 1(b) 
are normal crossing (the second subfigure shows the same configuration from 
a different angle with the surfaces clipped). Finally, in Fig. 1(c) the varieties 
are not normal crossing, because at the origin there is no regular system of 
parameters which could define both in the sense of Definition 2. 

3 The Hypersurface Case 

We use the notation and the context from Sect. 2, where X \, . . . , X s are pairwise 
different hypersurfaces (reduced separated subschemes of codimension one) . The 
following proposition provides an easy decision method to test normal crossing 
at a point. 
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(b) x 2 + y 2 + z 2 - 1, x 2 - y + z 2 , x 3 + 
2 



(c) x 3 — 3 xy 2 — z, 
(x,z) 



Fig. 1. Subfigures (a),(c): not normal crossing, (b): normal crossing 

Proposition 1. The hypersurfaces X\,...,X a C W, defined by fi,...,f s £ 
Ow respectively, are normal crossing at a point a £ W iff a £ Sing (A*) for all i 
and the {(dfi) (a) \ fi(a) = 0} are linearly independent. 

Proof. First let us note that f^w/k is a fc-vector space. Since the hypersurfaces 
that do not contain a are defined by the empty parameter system, we just con- 
sider those i for which a £ Xi, so that fiOw.a £ m a . If all these Xi are non- 
singular at a, the “gradients” ( dfi) (a ) are nonzero, and if they are also linearly 
independent, then fiOw,a form a regular sequence in Ow,a- This regular sequence 
is then extendible to a regular system of parameters. 

The “only if” part is immediate (see Definition 2 and Remark 2(ii)). 

The content of this proposition is well known, and the arising algorithm is the 
following. Let d Pl , . . . , d Pn be the derivations on W coming from a regular system 
of parameters. For a given a £ W, collect the /,; that vanish at a; without loss of 
generality we can assume that the relevant indices are 1 (< s). Compute 
defining equations for Sing(JC) as in Sect. 2 and check whether a £ Sing(Xj). 
If a is a singular point of any of the hypersurfaces, they are not normal crossing. 
Then compute (3w(/i, • ••,/;)) (a) and compute its rank. If it is less than l, the 
hypersurfaces are not normal crossing at a, otherwise they are. 

The next step is to generalize this method so that we can perform the test 
over closed subsets of W in one stroke of computations. It is quite obvious 
what we have to do. Instead of evaluation we check for a selected subset of 
hypersurfaces, say Xi, . . . ,Xi, whether there is a point a in their intersection 
where, after evaluation at a, we would get dfi, . . . ,dfi linearly dependent. 

We should remark in advance that this computation has to be performed 
only when l < n, since otherwise the arising matrix cannot be of rank l. Thus, 
when l > n the test reduces to check whether Xi fl • • • D Xi is empty, which 
corresponds with the fact that in an n dimensional ambient variety at most n 
hypersurfaces, intersecting at a, can have normal crossing at a. 

Assuming l < n, we have that 3iw.i(.f\ , ■••,/;) defines the points where the 
vector field given by dfi , ... ,dfi is linearly dependent (in the sense explained 
above). Nevertheless, only those points are relevant for us which are contained 
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in Xi fl • • • fl Xf, therefore the points in the intersection of the selected hyper- 
surfaces at which they are not normal crossing are defined by 

(/l) ■■■ifl) + 3W,j( fu • • • > fl)- 

We can immediately note that when the selected subset of hypersurfaces is 
a singleton, the above computation provides the singular locus of the hypersur- 
face. 

We can test now whether a given finite set of hypersurfaces in W are normal 
crossing by going through all the subsets, checking the selected hypersurfaces 
along their intersection. We use the method of Grobner bases (see [9, 10]) to test 
whether a given set of polynomials has a common zero. The algorithm isNCH 
summarizes our observations. 



Algorithm isNCH ; 

Input: W a nonsingular affine variety, (/i, ■ • • , /s) a list of defining 
equations of hypersurfaces in W. 

Output: true or false, indicating the normal crossing property of the input 
hypersurfaces in W. 
isNCH {W ,(/i , . . . , 

n := dim(W); R/J := Ow\ O := tdeg(Generators(7?)); 

for l = 1 to s do 

for each C {1, . . . , s} do 

if l < n then 

G := gbasis((/i 1 3W,i(/*i i - ■ • i fn) + J, 0)\ 
else G := gbasis((/ il , . + J, O); 

if NormalForm(l, G . ()) =/-- 0 then return false; 
return true; 



Remark 3. (i) The termination of the algorithm is easy to see. (ii) The correct- 
ness is supported by Proposition 1, which implies the correctness of the elemen- 
tary test for a selected subset of the input hypersurfaces, and the fact that the 
algorithm performs an exhaustive check, (iii) The algorithm is of exponential 
complexity in the number of input hypersurfaces, but please note that Grobner 
basis computation has double exponential worst case complexity in the number 
of indeterminates. 

4 The General Case 

We use the notation and the context from Sect. 2, and now X±, . . . ,X S denote 
pairwise different varieties (integral separated subschemes of W). We can easily 
check whether the varieties are nonsingular (see Sect. 2), thus till the formulation 
of the final algorithm we will just assume this property for all of them. 

For varieties of higher codimension, it is not enough to check linear depen- 
dence of vectors in their normal bundles. At a point a £ W, subvarieties with 
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normal spaces generated by linearly dependent vectors can either be normal 
crossing (think of a line embedded in a plane both embedded in 3-space) or not 
(if, for instance, they happen to be hypersurfaces). 

Therefore our strategy is to try to construct regular parameter systems which 
will serve as witnesses for Definition 2. We do it in two steps: first we achieve 
that each variety is defined locally by regular parameters (this is what we call 
a local regular representation, see Subsect. 4.1), and then we try to unify these 
parameter systems (to find local witnesses, see Subsect. 4.2). We must, of course, 
take care that, as in algorithm isNCH , we test all varieties against each other at 
all of their intersection points in W, see Subsect. 4.3. 

4.1 Local Regular Representations 

A nonsingular subvariety of W is a local complete intersection (see e.g. Exam- 
ple 8.22.1 on page 185 in [11]). The next definition just collects the data that 
represent such a variety by codimension many regular parameters. 

Definition 3. Let X C W be a nonsingular variety of codimension r. A triple 
( U,p , t) is a local regular representation of X (over U) iff U C W, p is a regular 
system of parameters over U and t £ {0, 1}”, such that X C\ U is defined by 
{tiPi | i = 1, . . . , n}. A list of local regular representations X — ((Ui,pyyU) \ i = 
1, . . . , s), is a regular representation of X iff X C (J? =1 Ui- 

Concerning the notation, elements of a regular system of parameters p are 
denoted by indexing: pi , while the zth parameter system in an enumeration is 
denoted by the parenthesized index: pyy The jth element of p(q is denoted 
by P(i)jy 

In this subsection we present an algorithm that computes a regular represen- 
tation for a given nonsingular variety. The algorithm is extracted from the imple- 
mentation of Villamayor’s algorithm for resolution of singularities by J. Schicho 
and the author (see [7]). To make the details of the computations readily avail- 
able first we briefly discuss two subalgorithms that will be used later. For the 
correctness proofs we refer to [7]. 

We do not need the matrix P of partial derivatives (see Sect. 2) in the local 
regular representations, but as partial derivation is ubiquitous in the algorithms, 
and we will have to change the parameter systems frequently, we present the 
transformation rules for P. If the reader is not interested in these details, they 
can be safely omitted. 

Given ([/, p, P), the first subalgorithm, Cover , passes to open subsets of U 
which are defined by its second parameter (/i, . . . , f s ) £ Ofj. The ith subset is 
isomorphic to U \ ZeroSet (/,;), which is described by introducing a new variable 
with a relation that makes it to be the inverse of fi. The p component does not 
change, while the P matrix needs to be updated for each subset, having a new 
column that contains the partial derivatives of the new variable with respect 
to p. The output subsets form a covering of U iff ZeroSet (/i, . . . , / s ) = 0. 

The subalgorithm Exchange replaces an element of the regular system of pa- 
rameters of (U,p, P) with a new function / £ Ojj, and it updates P accordingly. 
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Please note that the input specifications of the Exchange algorithm require that 
d Pi f is invertible in Ou- 



Algorithm Cover; 

Input: ( U,p , P), where U C W open, and p is a regular system of 
parameters over U with derivation matrix P ; and fi, . . . , f s £ Ou- 
Output: ((Ui,p;i), Pi) | i = 1, . . . , s) such that Ui = U \ ZeroSet (fi) and p 
is a regular system of parameters over Ui with derivation matrix Pi. 
Cover((U,p,P),(fi,-..,fs)) 

R/J := k[x i, . . . ,xi]/J := Ou ; 

for i = 1 to s do 

let Ui be defined by Ou t := R[xi+i]/{J + ( xi+ifi - 1)); 
let Pi £ Af £ xi+1 with ( Pi)mj ■= Pmj, for m = 1, . . . ,n, j = 1, . . . , l and 
(Pi)m,i + 1 := -x’i +1 d Pm fi, for m = 1, . . . ,n; 
return {{Ui,p, Pi) | i = 1, . . . , s); 



Algorithm Exchange; 

Input: ( U,p , P), where U C W open, and p is a regular system of 
parameters over U with derivation matrix P; i an index to p\ 
and /, q £ Ou such that q = {d Pi f )~ 1 . 

Output: (U,p ' , P') where pi is replaced by / in p' so that p' is a 
regular system of parameters over U with derivation matrix P' . 
Exchange{{U,p,P), i, /, q) 

{pi, ■ ■ ■ , Pn) := p; l := NumberOfColumns(P); 
p' := (pi, . . . , Pi-i, f ,Pi+i, ■ ■ ■ ,Pn)', let P' £ Ml xl ; 
for j = 1 to l do 

Pij := QPij) 

for m = 1 to n do if m ^ i then P) n j := Pmj — qPijd Pm f ; 
return (U,p ' , P')\ 



Let now X C W be a nonsingular subvariety of codimension r, defined by its 
vanishing ideal / = (fi, , fi)Ow- Let W be equipped with a regular system 
of parameters p with derivation matrix P. To construct a local regular represen- 
tation of X we intend to find elements in I that we can exchange by elements 
of p. We can have the following situations: 

1. There is a pi € I, so we just tag p 

2. We are not in case 1, but there is an / £ / and a p; € p such that d Pi f is 
invertible in Ow- I n this case we call Exchange to replace p; by / and we 
tag the new parameter /. 

3. We are neither in case 1 nor in case 2, but since X is nonsingular, there can 
be no point of X where all the partial derivatives of the defining equations 
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vanish. Then we can call Cover for (W,p, P) by ( d Pi fj \ i = 1 = 

1 , ,1). The subsets generated this way cover X , and over each of them we 
can call Exchange (since we made a partial derivative invertible), tagging 
the new parameters. 

This strategy is applied recursively over the subsets created by Cover , always 
considering only untagged parameters. The complete description can be found 
in the algorithm Regular Representation. 



Algorithm Regular Representation; 

Input: ( U,p,t ), P, where U C W open, p is a regular system of 

parameters over U with derivation matrix P, t £ {0, l} n tags for p ; 
and X a nonsingular subvariety of U. 

Output: A regular representation of X over U with a derivation matrix 
for each regular parameter systems. 

Regular Representation{{U , p,t), P, X) 

let I = {tiPi | i = 1, . . . , n) + (/i ,...,/;) define X such that 

0 mod {Upi | i = 1, . . . , n), j = 1, . . . ,Z; r := codimw(A); s:=^' I =1 t i ; 
/ / if codimension many regular parameters are tagged, we are done 
if s = r then return ((U,p,t), P); 

/ / case 1: look for a taggable p element 
for i = 1 to n do if ti = 0 and pi € I then 

ti := 1; return Regular Representation((U, p,t), P, A); 

// case 2: look for exchangeable elements over U 
for i = 1 to n do if ti = 0 then for j = 1 to l do 
if (dpi/j) -1 £ Ou then 

( U,p,P ) := Exchange ((U,p,P), i, fj, (dpjj) -1 ); U := 1; 
return Regular Representation((U, p,t ) , P, A'); 

// case 3: in general we have to pass to open subsets 
{{Ui,P(i),Pi)\i = 1 ,...,(n-s)l) := Cover((U,p, P), 

(9 Pi fj\j = l,...,l,ti = 0,i = l,.. . ,n)); 

LRR := (); N := 0; 

for i = 1 to n do if ti = 0 then for j = 1 to l do 

N := N + l;t' := t; t ' := 1; 

if d Pi fj ^ 0 then //we skip empty subsets 

( Un,P(n),Pn ) := Exchange((U N ,P( N ), Pn), i, fj, (dpjj)^ 1 ); 
append Regular Representation((UN ,P(N),t'), Pn, A) to LRR; 

return LRR ; 



Theorem 1. Let X C W be a nonsingular subvariety of codimension r , and letp 
be a regular system of parameters over W with derivation matrix P. The proce- 
dure call RegularRepresentation((lT,p, (0 | i = 1, . . . n)), P, X) computes a regu- 
lar representation for X in finitely many steps. 

Proof. The computation branches each time the Cover procedure is called. Thus 
we have to show that the algorithm computes a local regular representation 
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for X over an arbitrarily chosen path of this computation tree and that the 
open subsets at the leafs cover X. We also note here that we will not count the 
last recursive calls in which the “s = r” condition is satisfied and the procedure 
returns immediately. To emphasize this we will talk only about relevant recursive 
calls, in which real computations are done. 

We proceed by induction on r (< n). In the case when X is a hypersurface 
defined by f £ Ow, n o relevant recursive call happens. The first two cases 
(numbered by 1 and 2 in our previous enumeration) are straightforward. In 
case 3, we have 1 S (/) + ( d Pi f \ i = 1, . . . ,n) because X is nonsingular, thus the 
open subsets cover X and we are in case 2 over them. 

Let us assume that the algorithm fulfills the claims for codimension r — 1 
and let X be of codimension r (> 1). First we observe that the relevant part of 
the recursion is exactly r-fold over the branches where U does not get disjoint 
from X . This is because the algorithm tags one regular parameter in one recursive 
call and a subvariety defined by i regular parameters is of codimension i. 

In the first call of the algorithm for X, if case 1 or case 2 applies, we imme- 
diately reduce the problem to solving one of codimension r — 1 over the closed 
Y C W defined by the tagged parameter. In case 3, considering X as a subvari- 
ety of the nonsingular subset of W defined by the tagged parameters, we have 
1 € (/i, . . . , fi) + (d Pi fj \j = l,...,l,ti = 0,i = l,...,n) again because this ideal 
defines a subset of Sing(X) (which is defined via the Jacobian ideal) and X is 
nonsingular. Thus the union of the computed subsets contains X and over the 
nonempty ones case 2 applies. Therefore over a subset U from these last ones, 
we again reduced the problem to one of codimension r — 1 over the closed Y C U 
defined by the tagged parameter. 

This reduction is possible because being nonsingular is an intrinsic property 
of X and because X is a local complete intersection in any nonsingular variety 
containing it. Please note that Y is nonsingular in U and that the untagged 
elements of p form a regular system of parameters over Y . 

Finally, by the induction hypothesis, the algorithm computes a regular repre- 
sentation for X over Y which extends to a regular representation of X (possibly 
via the open cover constructed in case 3) with the tagged parameter. 

Remark C (i) A rough worst case complexity bound for the algorithm is ( nl) r 
subsets in the final regular representation, where l is the maximal cardinality 
of the generating sets of the defining ideals used for X throughout the algo- 
rithm, and X is of codimension r in the n-dimensional W. (ii) The correctness 
proof does not talk about P, since this datum is not a member of a regular 
representation. However the correctness of these matrices, attached to the local 
regular representations of the output, can be easily seen via the correctness of 
the subalgorithms Cover and Exchange. 

4.2 Common Regular Parameter Systems 

We consider the original scenario, where we have Xi , . . . , X s , subvarieties of W. 
Now that we can compute regular representations for each of them, our next 
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goal is to test local “unifyability.” More precisely, if a selected subset of the given 
varieties, say X \, . . . , Xi, has a nonempty intersection S , we want to determine 
whether for each point of S there exists a neighborhood U together with a regular 
system of parameters such that each selected variety is defined by a subset of 
this parameter system over U (see Definition 2). 

In this section we use the notational conventions introduced in Definition 3, 
and in the paragraph after it. 

Definition 4. Let P(i), • • • iP( s ) be regular parameter systems, (i) We call an 
element f £ Ojj locally exchangeable at index j for P(i), ■ ■ ■ ,P( S ) iff dp^jf 0 
for all i = 1, . . . , s. We call f locally exchangeable over some open subset V C U 
(or at a point a £U) ifZeroSet(dp^ ):j f)nV = 0 (or (<9 P(i)j . f)(a) yf 0 respectively) 
for alii = 1 (ii) Arrangements of Pm, • • • ,P( S ) are elements of the set 
of all possible permutations of the indices of the parameters in the systems, i.e. 

{ (.^lP(l) > • • ■ j &sP(s)) | CTi , . . . , (J s £ S'n} • 

Remark 5. (i) An element of Ow which is locally exchangeable at a £ W neces- 
sarily defines a nonsingular hypersurface at a. On the other hand, if an element 
defines a singular hypersurface at a, it cannot be exchangeable over any neighbor- 
hood of a at any index with respect to any parameter system, (ii) The property 
of being locally exchangeable is locally stable, i.e. an element with this property 
for a given list of regular parameter systems and index j at a point a £ W has 
this property at all points of some neighborhood of a (cf. Remark 1). 

Example 1. Let W = A |, p (!) = (xi, x 2 — x\+x\, x\+x 2 — xf), p( 2 ) = (xi, x 2 , x 3 ) 
and let / = x 2 . The partial derivatives of / with respect to p(i) are (3x^ — 1, 1,0), 
thus / is locally exchangeable for p^ at index 1 and 2. The corresponding open 
subsets over which the exchange is possible are W \ ZeroSet(3xf — 1) for index 1 
and W for index 2. 

The partial derivatives of / with respect to p( 2 ) are (0,1,0), thus / is ex- 
changeable for p( 2 ) only at index 2, actually we have P( 2 ) 2 = /. 

If we have to keep P(i) 2 in p(i) because it is used in a local regular repre- 
sentation, and we need to exchange / for p^ at index 1, we have to rearrange 
either pm or p( 2 )- Say, we keep p^) fixed, and take <r 2 to map (1, 2, 3) i— > (2, 1, 3), 
getting <r 2 p( 2 ) = (x 2 ,xi,X 3 ). In this case we have that / is exchangeable for the 
arrangement p(i), <J 2 P( 2 ) at index 1. 

The next lemma allows us to restrict the search for exchangeable elements 
in an ideal to a finite generating set. 

Lemma 1. Let P(i), ■ • • iP(s) be regular parameter systems over W. The vanish- 
ing ideal I C Ow of an algebraic set X C W contains an element which is locally 
exchangeable at index j for Pm, ■ ■ • ,P(s) over some neighborhood of a £ X iff 
any finite generating set of I does. 

Proof. The “if” part is immediate. For the “only if” part, let us assume that 
there exists / £ I locally exchangeable at index j for P(i), ■ ■ ■ ,P( S ), while there 
exists a finite generating set of I in which no element has this property. Now 
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observe that if all the elements from the finite generating set are not locally 
exchangeable then neither are their Ow~ linear combinations. We just have to 
apply the Leibniz rule and the fact that both the generators and their partial 
derivatives, for some of the parameter systems with respect to the jth parameter, 
vanish at a. Therefore, as we must be able to write / as a linear combination of 
the generators, / also cannot be locally exchangeable. 

The next lemma provides the core idea of the algorithm that tries to find 
a common regular representation for varieties from their given regular repre- 
sentations. It shows that locally at a point normal crossing of the varieties is 
equivalent with the existence of suitable arrangements of the parameter systems 
such that at each index there is an exchangeable element in the intersection of 
the ideals of the relevant varieties. Relevance just means that the corresponding 
local regular representation has a (1) tag at the index in question. 

Lemma 2. Let (Ui,p(i),ti) be a local regular representation of Xi C W for 
i = 1, . . . , s, let L C Ojji define Xi over Ui and let a £ fji=i(^i LI Xf). There 
exists an arrangement of the the parameter systems , and of the values in the 
tags accordingly, such that for each index j there exists pj £ flt ( } =i (th e 

empty intersection defines Ow ) with the property that pj is locally exchangeable 
at index j for all the p ^ over some neighborhood of a and pj — pj (a) fL (pi — 
Pi (a), . . . ,pj - 1 — pj-\(a))Ow,a when j > 1 iff the Xi are normal crossing at a. 

Proof. Let us assume that the X, are normal crossing at a, so that we have a reg- 
ular system of parameters p such that each X, can be defined by p-parameters 
locally at a. For each i we can choose arrangements of pu\ for which each pj 
is locally exchangeable at index j for po) with the following properties: (i) pj 
becomes a parameter over some neighborhood of a. This can be achieved by en- 
suring that (dp ({ . .pj)(a) ^ 0, and the existence of a suitable arrangement follows 
from the fact that p and p(q are regular parameter systems over neighborhoods 
of a. (ii) When j > 1, pj - pj(a ) fL (pi -pi(a), . . . ,pj - 1 - pj-i(a))O w ,a, which 
follows immediately from the fact that p is a regular system of parameters over 
some neighborhood of a. 

Now since p is a witness of normal crossing for the varieties at a, Xi can be 
defined by codimension many p-elements at a. Let these elements be identified by 
the tag Ti. But p(q is also a regular parameter system with the same property, so 
there must be an arrangement of p(q which fulfills requirement (i) and (ii) above 
and for which T; = L, otherwise the tagged parameters could not define the same 
variety at a. Choosing such an arrangement for each p ^ we get pj £ (~| t =1 Ii. 

For the “only if” part, let us assume that the parameter systems are rear- 
ranged so that there exists pj £ p| t , =1 h with all the required properties for 
j — 1 , ,n. Then p\, ... ,p n forms a regular system of parameters over some 
neighborhood of a such that Xi is defined locally at a by the parameters tagged 
by ti . That is, the Xi are normal crossing at a. 

Next we informally present the algorithm that tries to prove normal cross- 
ing of given varieties along their intersection by constructing witness regular 
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parameter systems. For given local regular representations (Ui,p(i),U) of the va- 
rieties Xi, i = 1, . . . , s, let U = P)®_i Ui ; we consider all the regular parameter 
systems and varieties restricted to U. Let define Xi over U and let us also 
assume that S = U fl fT- =1 Xi is nonempty (otherwise we have nothing to do). 

We try to unify the using the observations of Lemma 2, searching through 
all arrangements of the parameter systems. Our first goal is to compute local 
regular representations in which the regular parameter systems are the same for 
each X, . Thus we start from a given list of local regular representations (one for 
each Xi) and search through arrangements of their parameter systems to find 
ones that allow exchanges of elements from the corresponding ideal intersections. 

Algorithm CommonRR; 

Input: A list of regular representations {(Uij,puj\, tij) \ j = 1, . . . , rrii) of 
varieties Xi, i = 1 , . . . , s, with derivation matrices Pij. 

Output: true, and ((({7L,pL tC) | i = 1 , . . . , s) \ j = 1, . . . , m) a list of 
branches giving rise to regular representations of the Xi over 
U' = au U ij A S — -Xi, where for each fixed j, P(%j ; is the same 
for all i. Or false if the varieties are not normal crossing over W, and 
then the list of regular representations is such that S\U' are the 
points where the varieties are not normal crossing. 

CommonRR(((((Uij ,p( i;j ),tij) , P^) \ j = 1, . . . , m/) | i = 1, . . . , s)) 

LW := ();$:= flLi 

if 5 = 0 then return (true, ()); 

for (ji,...,j s ) = (1, — , 1) to (mi,...,m s ) do 

for each arrangement of the parameters in PnjPj, ■ ■ ■ iP(sj s ) do 

L := (JJSJJiji , P(iji) , tiji ) , Fiji ) l* = 1 , - - ■ , s)) ; // start with a single branch 
for j = 1 to n do N :=();// N collects branches created at index j 
for each (((f7i,P(*),t*),Pi) | i = 1, . . . , s) in L do T := (h \ t Wj = 1); 
u := H- = 1 Ui-, S' = UnS] let h C O v define W; F := (); 
if #T = 0 then F := (p {h)j \ d Pmi p Wj j- 0, l = l,...,s); 
elif #T = 1 and d mj p( T[1])j ± 0, l = l,...,s then F := (p( T [i])j); 

else </i, D t=i It[X\O v \ F := (f h \ S' % Sing(ZeroSet(/ ft )), 

fh 7^ 0, l = 1 , . . . , s, h= 1 , . . . , r); 

if = 0 then next; // this branch is dropped and we take the next 
/ / each element of F will give rise to a separate branch 

for h = 1 to #F do 

//we get rid of extraneous components and prepare exchange 
for i = 1 to s do (/i, . . . ,/ n ) := {{t{i)jF[h]) + ft^iP(i)i 1 1 ± j)) ■ h\ 
((Vf,p{ i) ,Pl)\l = l,...,r i + 1) := 

Cover((Ui,p {i) ,Pi), ( F[h\ ) + (f m d P(i)J F[l\ \ m= 1 ,...,n)f, 
for (h, . . . , l s ) = (2 , . . . , 2 ) to (n + 1 , . . . , r s + 1) do 
/ / exchange and create new branches 
for i = 1 to s do (V /‘ , p ^ ,P/) := 

Exchange ((Vf, p 1 ^ ,P^), j, F[h ], (d^F^])" 1 ); 
append (((Vf,p l fo, U), P}*) \ i = 1, . . . , s) to N; 
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L := N- 

append (((£/*, p w , U) \ i = 1,. . . , s) | S' D f|i=i u i ± 
muPM,U),Pi) \i = e L) to LW\ 

if S C LW\j}[i][l] then return (true, LW)' 

return (false, LW)', 

Such a list of local regular representations that we handle together in the hope 
that we will be able to reach the first goal (via several parameter exchanges and 
maybe by passing to open subsets) we call a branch. 

Definition 5. Using the notation of the previous three paragraphs , we define a 
branch to be a list {{Ui,pu),ti) \ i = 1, . . . , s), where (Ui,p^\,ti) is a local regular 
representation of Xi . We call a branch a local witness of the normal crossing of 
the Xi along their intersection in U iff all the p ^ in it are the same over U . 

Our second goal is to find sufficiently many local witnesses such that the in- 
tersection for all i of the open subsets over which the local regular representations 
give rise to regular representations for the Xi, contain S. 

We organize the search by going through all arrangements of the parameter 
systems, treating the indices incrementally. We may have to pass to open subsets 
in order to perform the exchange operations, which results in a tree of branches 
as the search progresses. We always consider the leaves of this tree. If on a certain 
branch we find no locally exchangeable elements at a certain index, the branch 
is dropped. Completed branches, on the other hand, provide local witnesses. 

The complete description of unifying local regular representations can be 
found in algorithm CommonRR. This time, in order to keep the presentation 
simple, we do not expand purely technical subalgorithms like “go through each 
arrangements of the parameter systems.” As usual the instruction “next” starts 
immediately the next execution of the innermost loop. 

Theorem 2. Given a list {{{{Uij,puj),tij), Pij) \ j = l,...,mj)|* = l,...,s) 
of regular representations for the Xi. If the represented varieties are normal 
crossing along their intersection, the procedure call CommonRR returns true 
and a list of local witnesses as given in its output specifications in finitely many 
steps. Otherwise it returns false, and the local regular representations cover some 
neighborhood of the points of the intersection of the Xi at which they are normal 
crossing but contain none of the points where they are not normal crossing. 

Proof. The notation comes from the algorithm CommonRR. The termination of 
the procedure follows from the fact that each loop runs only finitely many times 
and in the innermost loop only finitely many branches are created at each run. 

Assume that the input describes varieties which are normal crossing along 
their intersection. The existence of suitable arrangements of the parameter sys- 
tems to construct local witnesses of normal crossing in some neighborhood of 
a point a € S is ensured by Lemma 2. The construction, i.e. the local regular 
representation, is locally stable, thus, because of noetherianity, there is always 
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a finite covering of (some neighborhood of) S so that one can choose the same 
witness at each point in an open subset of this covering. 

Since the algorithm examines all possible arrangement of the parameter sys- 
tems starting with each possible initial branch and considers each locally ex- 
changeable element from finite generating sets of the corresponding intersec- 
tions of ideals at each index, which is sufficient by Lemma 1, moreover because 
it takes the largest possible open subsets over which the parameter exchange can 
be performed, the resulting local regular representations cover S. 

There are two more details to be settled here. The first is that since we have 
to exchange elements into the parameter system of a local regular representation 
of a variety that come from the intersection of its ideal with other ideals, after 
exchanging the parameter, extraneous components may appear. We eliminate 
them by passing to open subsets which do not contain them and whose union 
contains all points in their complement. 

The extraneous components cannot meet S because by definition the new 
parameter (which introduced them) has normal crossing with all the other pa- 
rameters in the new system, at each point of the intersection of S and the open 
subset over which the exchange is possible. If an extraneous component passed 
through one of these points, the variety would become singular there and this 
would contradict the fact that the new parameter system provides a local regular 
representation for it at the point in question. 

The other issue is to make sure that for each j the new parameter pj ful- 
fills pj —pj(a) $ (pi —pi(a), . . . ,Pj~ 1 — pj-i(a))Ow,a for all a in the open subset 
over which the new parameter system will be considered. This is ensured by the 
fact that at index j we test local exchangeability with respect to the parameter 
systems in which the first j — 1 entries are already the new ones. That is, we con- 
struct the regular parameter system of a local witness (which is the same for all 
the local regular representations in it) in a way that we have regular parameter 
systems at each intermediate step of the process. 

Assume now that the varieties are not normal crossing over S C S, and let us 
take an a £ S. By Lemma 2 at each arrangement of the parameter systems there 
is an index j at which we cannot find a locally exchangeable element in the inter- 
section of the ideals of the varieties identified by the tags at index j. Therefore, 
by Lemma 1, no branch containing a can complete and enter in LW. On the 
other hand, the varieties are normal crossing over W \ S thus the construction 
applies for this subset and the output verifies the normal crossing property over 
some neighborhood of S\S. 

Remark 6. The worst case complexity of the algorithm is clearly super expo- 
nential. This causes a problem from the efficiency point of view mostly when 
the varieties are not normal crossing, since this fact is found out only after the 
search is completed. To ease this problem, if one does not need the regular rep- 
resentations over S\S (defined in the proof of Theorem 2) , in case the varieties 
are not normal crossing, the algorithm can be interrupted as soon as an ideal 
intersection is found in which all the (finitely many) generators define singular 
hypersurfaces at some point of the intersection of the varieties (see Remark 5(i)). 
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There is another technique with which we can reduce the number of local 
witnesses in the output and which deserves a few more words. It is referred by the 
name “focus” strategy in [12, 13]. Briefly, it requires the addition of a new data 
entry to local regular representations, let us call it F, which is a (possibly empty) 
list of polynomials. Consider a regular representation X of some nonsingular X. 
The role of F is to define the subset of C7, of a local regular representation 
in X, which is not covered by the Uj of X with j < i. Initially F := () and in 
the exchange operations it does not change its value. In a cover operation with 
input /i, . . . , fi, when i > 1, for the ith new subset we append /i, . . . , /)_ i to F 
(which is the “focus” for the input subset), to obtain the new “focus” for the ith 
output subset. 

We benefit from the data of F in the computations in that we know, when 
ZeroSet(P) = 0 over the underlying subset Ui , every point of Ui is already 
contained in subsets with smaller indices, thus Ui is redundant. Or, since for us 
only S , the intersection of the varieties is interesting, we can declare Ui redundant 
even if ZeroSet(P) D S = 0. If in a branch we find that all the Ui are redundant, 
the whole branch can be dropped without loss of data that is relevant for the 
computation. 

4.3 Normal Crossing Test for Varieties 

The test of the normal crossing property for given varieties X\, . . . , X s C W is 
described by the algorithm isNC. It uses the algorithm RegularRepresentation 
to compute initial representations for the varieties. Then for each subset of the 
set of input varieties, the corresponding regular representations are checked by 
the algorithm CommonRR along their intersection. 



Algorithm isNC ; 

Input: W a nonsingular affine variety, and varieties Xi, . . . , X 3 C W. 
Output: true or false, indicating the normal crossing property of the input 
varieties in W, and a list of branches giving rise to regular 
representations of the X, over some neighborhood of the points 
where the X, are normal crossing. 
isNCH(W , (X 1 ,...,X S )) 

let p be a regular system of parameters over W with derivation matrix P; 
n := dim(tr); NC := true; 

for l = 1 to s do 

Ri := Regular Representation^, p, (0 | i — 1, . . . , n), P, V); 

for l = 1 to s do 

for each C {1, . . . , s} do 

(N, L) : = CommonRR ( (Pi 1 , . . . , Pi, ) ) ; append L to LW ; 
if not N then NC := false; 
return (NC, LW)', 
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Theorem 3. Given Xi,...,X s C W, calling isNC (W,(X i, . . . ,X S )) deter- 
mines in finitely many steps whether the varieties are normal crossing over W 
and returns a list of local witnesses which gives rise to regular representations 
of the varieties over some subset of W which does not contain the points where 
the Xi are not normal crossing. 

Proof. The claim follows from the termination and the correctness of the subal- 
gorithms and from the fact that the algorithm performs an exhaustive search. 

5 An Extended Example 

The algorithm in Sect. 3 is straightforward; therefore in this section we present 
an example only for the general case, discussed in Sect. 4. 

Let k := Q, W := A |, O w ■= k[xi,x 2 ,x 3 \, p := (xi, x 2 , x 3 ), P := 1 3 . 
Let Xi be defined by Ji := (x 2 , (x 4 — l) 2 + x 3 — 1) and X 2 be defined by I 2 := 
(x 3 , (xi + l) 2 + x\ - 1) (see Fig. 2(a)). 

In the first phase, we compute local regular representations for the varieties, 
e.g. for X\ we call RegularRepresentation((W,p, (0, 0, 0)), P, Xi). The algorithm 
finds that x 2 £ I± is a regular parameter (case 1); therefore it tags it, and calls 
itself with (W,p, (0, 1, 0)), P, X 1 . Naturally h = ( p 2 ) + (/i) = ( x 2 ) + ((xi - 
l) 2 + x| — 1) where {x\ — l) 2 + x\ — 1 ^ 0 mod x 2 , and we are in case 3 since 
neither d Pl f\ = 2xi — 2 nor d P3 f\ = 2x 3 is invertible over W. 

We pass to open subsets by Cover((W,p, P), {d Pl fi,d P3 fi)), Obtaining: 

U\\ Ojj 11 = k[x\, x 2 , x 3 , X4]/(x4(2xi — 2) — 1), over which we can call 
Exchange((Un , p^i) , Pn), 1, (xi - l) 2 + x\ - l,x 4 ), 

U 12 Ou 12 = k[x\, x 2 , x 3 , x 4 ]/( 2 x 4 X 3 — 1), over which we can call 
Exchange((Ui 2 ,p( 12 ),Pi 2 ),3,(xi - l) 2 + x\ - l,x 4 ). 




( a ) (b) (c) 



Fig. 2. (a): Xx,X 2 , (b),(c): the relevant subsets of the local regular represen- 
tation of X\ 
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Over U\\ and Ui 2 the situation is analogous (see Figs. 2(b)(c)), so let us take, 
say, C/11, where after the exchange operation we get (f7n,p(n) = ((*1 — l) 2 + 
xl - 1, x 2 , x 3 ), (1, 1, 0)) with P u ■■ 

/ X4 0 0 — 2^4 \ 

0 10 0. 

\-2£ 3 £ 4 0 1 Ax 3 x%) 

In the following we do not present the derivation matrices since they become 
more cumbersome. 

Over U\\ the first parameter P(n)i gets tagged while over C/ 12 it is the 
third, P(i2)3- For these subsets the next recursive call will find s = r, i.e. the 
number of tagged parameters is the same as the codimension of X \ , so it returns 
immediately. 

The computed regular representation X\ of Xi consists of two local regular 
representations: 

- (C/n, (On - l) 2 + xl - 1, *2, *3), (1, 1, 0)), 

- (C12, (x!,x 2 , (On - l) 2 + x\ - 1), (0, 1, 1)). 

Analogously, for X 2 we get X 2 : 

- (C/21, (On + 1) 2 + x%-i, x 2 , x 3 ), (1, 0, 1)), 

- (C/22, On, (n + 1) 2 + x \ - 1, x 3 ), (0, 1, 1)). 

Then, in the second phase of the algorithm, for each nonempty subset of 
{Xi, X 2 } we call CommonRR. The first is {X 3 }, in which the intersection S is Xi 
itself, and the outermost loop simply runs through the local regular representa- 
tions in X\. Each of them will result an execution without branching, because 
the “common” parameters are the ones in the regular parameter systems no 
matter which arrangement we take. When the second local regular representa- 
tion is processed for the first time the condition that S is covered by the [/, of 
the local regular representations in LW is satisfied and the algorithm returns 
true. The procedure is analogous for {X 2 }. 

When {Xi, X 2 j is passed to CommonRR, nothing interesting happens when 
one of the chosen local regular representations is over Ui 2 or U 22 . In these cases 
one (or both) of the varieties is outside U , the open subset under consideration. 

So let us take (Uu, ((£1 — l) 2 + £§ — 1, x 2 , x 3 ), (1, 1, 0)) for Xi and (C/21, ((£1 + 
l) 2 T x 2 — 1 , x 2 , x 3 ), ( 1 , 0 , 1 )) for X 2 . Their intersection U can be defined by Ojj = 
k[xi,x 2 ,x 3 ,X4,x 3 \/ {xi(2xi — 2 ) — 1 , £5(22:1 + 2 ) — 1 ), and the inclusion map 
U Uu gives rise to the ring homomorphism Ou, — > Ojj, %i Xi, i = 1, 2, 3 , 4 , 
while for U 2 i it is Xi 1— > Xi, i = 1 , 2, 3 and £4 1— > £ 5 . 

For the arrangement that does not change the initial order of the parameters, 
we get that for index 1 there are two tagged parameters and that in the ideal 
intersection ( Ii n I 2 )Ou there is an element which is exchangeable for both 
parameter systems at this index: / = x\ + x%xi + £i£§ + 2£ 2 — Axi — 2x\ (see 
figure 3(a)). The redundant component appearing in Xi is a line parallel to 
the £3-axis standing at the point (—2,0,0), whose ideal over U is generated 
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(a) (b) (c) 



Fig. 3. (a): ZeroSet(/), (b),(c): extraneous components appearing in Xi,X 2 

respectively 



by x 2 ,xi + 2 (see Fig. 3(b)). Moreover the element we have to make invertible 
over Un is 9 P{11)1 / = (3xf + x\ + x 3 — 4)2:4, which is nonvanishing at S. 

For X 2 the redundant component is another line parallel to the a^-axis pass- 
ing through the point (2, 0, 0), whose ideal over U is generated by 2:3, X\ — 2 (see 
Fig. 3(c)). Over U21 the partial derivative which we have to make invertible is 
formally the same as the one over Uu (as / is symmetric in x 2 , X3 and swapping 
them in P(u) 1, P(2i)i induces 2:1 —x±, which is hidden in the partial derivative 

since it is quadratic in X\): <9 P(2)1 / = (32:^ + x\ + x\ — 4)2:4. 

In the cover operations we have to make the corresponding partial derivative 
and a generator of the extraneous component invertible simultaneously, thus we 
take their product. For instance Uu is covered along /, x 2 d P{11)1 f , (xi+2 )<9 P(11)1 /- 
The first subset will always be disjoint from S ; therefore it will be discarded in 
the next loop, where the subsets obtained from Uu and U 2 i are combined. In 
this example four branches are created and collected in N. 

After the exchange operations are done, in the new branches we have regular 
representations for X\ and X 2 in which the first regular parameters are the same. 
As the algorithm continues for index two and three, it needs to do nothing since 
there is only one parameter tagged at both indices which actually appears in 
both of the parameter systems. Thus the branches complete, and since S will 
be contained by the intersection of the union of the open subsets in the regular 
representations of X\ , X 2 given by these branches, the termination condition is 
fulfilled. 

We can also note that the branches over which x 2 S P(11)1 / or £3 <9 P(21)1 / were 
made invertible are again disjoint from S. Thus the only interesting branch is: 

- X i: x 2 , x 3 ), (1,1,0)), 

- A 2 : x 2 ,x 3 ), (1,0,1)), 

where O v = k[xi, x 2 , x 3 , 2:4, 2:5] 7(2:4 (22:1 - 2) - 1, 215(2:1 + 2)(<9 P(11)1 /) - 1) and 
O v = k[ 2:1, 2; 2 , 2:3, 2:4, 2:5]/ (2:4 (2xi + 2) - 1,2:5(211 - 2)(d P(21)1 f) - 1) (the partial 
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derivatives are naturally included in the superrings) . This branch is a witness of 
the normal crossing of X\ and X 2. 
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Abstract. In this paper, we present an algorithm to compute the pro- 
jection of a quasi variety over an algebraic closed field. Based on the 
algorithm, we give a method to prove geometric theorem mechanically, 
and the non-degenerate conditions that we get by the method are proved 
to be the ’’weakest”, i.e. the geometric theorem is true if and only if 
these non-degenerate conditions are satisfied. A method for automatic 
geometric formula deduction is also proposed based on the algorithm. 
The algorithm given in this paper has been implemented in computer 
algebra system Maple. 



1 Introduction 

Wu’s method [1], as well as others based on computing of Grobner basis, has 
been fruitfully applied to automatic theorem proving in elementary geometry. 
Using Wu’s method, more than 600 theorems have been proved by computer [2]. 
For both Wu’s method and Grobner basis method, often a geometric theorem 
is true only in a ’’generic” sense, that is, certain degenerate cases must be ruled 
out. 

As for other methods based on Grobner bases, various way to find non- 
degenerate conditions are reported. Kapur [3] describes a complete method of 
Grobner bases to prove geometry theorems, including how to obtain non- 
degenerate conditions. As Kapur claims that ’’conditions found using this ap- 
proach are often simple and weaker than the ones reported by using Wu’s method 
or reported by Kutzler & Stifter’s [4] and Chou & Schelter [5] based on the 
Grobner basis method.” 

F. Winkler [6] has given a method whereby to compute the simplest non- 
degenerate conditions. His criteria for simplest is ”of as low a degree as possible” 
or ’’involving only certain variables.” 

In many cases, the non-degenerate condition is too strong. The theorem is 
still true even the non-degenerate condition is not satisfied. 



F. Winkler (Ed.): ADG 2002, LNAI 2930, pp. 21-30, 2004. 
(c) Springer- Verlag Berlin Heidelberg 2004 
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In this paper, we present an algorithm of the projection of quasi variety 
which is based on Wu’s method and his theorems about the projection of quasi 
variety [7]. A new method is given for mechanical geometric theorem proving and 
formula deduction based on the projection of quasi varieties. During proving 
geometry theorem, the weakest non-degenerate conditions can be obtained by 
this method. Moreover, the variables occurring in such conditions are determined 
in advance. As for geometry formula deduction, we can derive all conclusions 
about certain variables. At last, we give two examples to show how the method 
works. 

2 Non-degenerate Condition and Projection of Quasi 
Variety 

In this paper, all polynomials are in polynomial ring K[x i, • • • x n ], where K is 
a computable field of character 0. E is an algebraic closed extension field of K. 
For a geometric theorem, after setting up an appropriate coordinate system, the 
corresponding geometric configuration of hypothesis can be expressed by a finite 
set of polynomial equations PS = 0 i.e. PS = {pi,--- ,p s },pi = 0, • • • ,p s = 

0 , pi £ K[xi,--- ,x n }. The geometric configuration of the conclusion can be 
expressed by a polynomial equation C = 0, i.e. c = 0,c £ K[x i,--- ,x n }. If 
PS = 0 => C = 0, then the theorem T = (PS, C ) is called to be universally true. 

1. e. \/x £ E n , p\(x) = 0, • • • ,p s (x ) = 0 => c(x) = 0. In most cases, a geometric 
theorem is not universally true. It is true only if the non-degenerate condition 
is satisfied, j ^ 0 is called the non-degenerate condition if 

(I) (Vx £ K n )(pi(x) = • • • = Ps (x) = 0 A g(x) 0 => c(x) = 0) 

(II) (3x £ K n )( Pl (x) = • • • = p s (x) = 0 A g(x) ^ 0) 

From the above definition, we can see that the non-degenerate condition is a suf- 
ficient condition for a geometric theorem to be true. 

Wu’s non-degenerate condition: In Wu’s method, the non-degenerate condition 
is defined as the product of the initials of the characteristic set not equal to zero. 
Suppose CS : , G m is the characteristic set of PS, J = ■ /* where 

is the initial of G^. If the pseudo-remainder of the conclusion polynomial G 
w.r.t the characteristic set CS is 0, i.e. there are non-negative integers Sj s.t. 
7* 1 • • • P^C = Q\C\ + • • • + QmCm + 0, then J ^ 0 is the non-degenerate con- 
dition for the theorem to be true. 

Kapur’s non-degenerate condition: Let G\, Gi be the Grobner basis of ideal (PS) 
and (PS U {Cz — 1}). if Gi ^ {1} and Gi ^ {1}, then the theorem is true 
when gi 0. If gi satisfy 

(a ) gi £ G 2 D K[x\ A g t £ (PS) 

(b) l ^ GrobnerBasis(P5' U {giZ — 1}) 
gi ^ 0 is the non-degenerate condition. 
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Winkler’s non-degenerate condition : Winkler has proved that all polynomials 
satisfying (I) constitute an ideal, and among the polynomials in the Grobner 
basis of the ideal which satisfies (II), there is a polynomial g which has the least 
leading term, 3 / 0 is the simplest non-degenerate condition. 

All the non-degenerate conditions mentioned above, are sufficient conditions 
for a geometric theorem to be true. 

Let PS be a polynomial set and Dbea single polynomial in K[x 1 , • • • , x n \ . E 
is an algebraic closed extension field of K as before. We define 

Zero(PS) = {e £ E n \\/P € PS 7 P(e) = 0} 

Zero{/D ) = {e £ E n \D(e) ± 0} 

Zero(PS/D) = Zero{PS ) f) Zero{/D ) 

Definition 1. ( Quasi Variety) For any finite number of polynomial sets PSi and 
polynomials G, in K[x 1 , ■ ■ • , x n ], the set 

U ZeroiPSi/Gi) 

i 

is called a quasi variety. 

Definition 2. (Projection) r To eliminate variables x m +\,-“ ,x n , the map pro- 
jection is 

Pr 0 jx m+ i,- - ,x n : E n — > E m 

which sends (oi,--- ,a n ) to (ai,--- ,a m ). If V is an affine variety in E n , 
P r °jx m + (V) may not be a variety in E m 7 but it is contained in a va- 
riety in E n . 

For a polynomial set PS and a polynomial D , we apply projection 
Projxm+ 1 ,- ,x„ to Zero(PS/ D). Then we have 

Proj Xm+1 ,..., Xn Zero(PS/D) = {e£E m \3a £ a) 

£ Zero(PS/D )} 

when m = 0, we define Proj Xlt ... iXn Zero(PS/ D) = true , if Zero(PS/ D) ^ 0; 
and false otherwise. 

The projection of a quasi variety in E n to E m is a quasi variety in E m . Please 
see[7] for the details. The projection of a quasi variety has been investigated by 
Wu [7], Wang [8] and Gao [10]. 

3 Algorithm to Compute Projection of Quasi Variety 

In Wu[7], a method for computing the projection of a quasi variety is given. 
Before giving a new algorithm, we will give several theorems which are needed 
for constructing the new algorithm. 
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Theorem 1. Let PS be a polynomial set and D a polynomial , there is an 
algorithm to decompose PS into a finite set of ascending set ASi s.t. 

Zero(PS/ D) — \^J Zero(ASi/ JiD) 

% 

where each ASi is an ascending set, Ji is the production of the initials of the 
polynomials in ASi. 

The proof and the algorithm can be found in Wu [1]. In [9], an improved algo- 
rithm has been given to decompose a polynomial set into a series of ascending 
sets. 

Lemma 1. If Zero(PSfD) = Zero(ASt/ JiD), then we have 

Proj Xm+lj ... , Xn Zero(PS/D ) = (J Proj Xm+l! ... , Xn Zero(ASi/ J t D) 

i 

Proof. It is enough to prove 

P ro jx m+1 ,... ,x n U; Zero(ASi/JiD) = U- Proj Xm+1} ... , Xn Zero(ASi/ JiD). 

Take any element a = (ai,--- ,a m ) £ E m from Proj Xm+1: . .. tXn (J- Zero(ASi/ 
JiD), then there exists a' = (a m+ i,--- ,a n ) £ E^ n ~ m ^ s.t. (ai,-- - , a ra )£lJi 
Zero(ASi/ JiD), then there exists an i s.t. (ai,--- ,a n ) £ Zero{ASi/ Jf). Accord- 
ing to the definition of projection, (ai, • • • , a m ) £ Proj Xm+lt ... tXn Zero(ASi/ Ji), 
it follows that a = (ai, • • • , a m ) £ UiProj Xm+li ... iXn Zero(ASi/ JiD). It shows 
Proj Xm+lr .., Xn IJi Zero(ASi/ JiD) C U i Proj Xm+u ... , Xn Zero(ASi/ JiD). 

The inclusion of reversal direction also can be proved by the same way. 

Lemma 2. AS = {Ai, ■ • ■ ^4 S } is an ascending set w.r.t. variable ordering (x\ < 
*2 < ••• < a : n ),AS' = {Ai,- - ,yl s _i}, J and J' are the products of the ini- 
tials of polynomials in AS and AS' respectively. I s is the initial of A s . d = 
degree(A s ,x n ), R = Prem(D d , A s ,x n ). 

1. if degree(A s ,x n ) ^ 0, degree(D,x n ) ^ 0 then 

Proj Xm+lt ..., Xn Zero(AS/ JD) = Proj Xm+u ... tXn Zero(AS' / J' (I S R)) 

2. if degree(A s ,x n ) ^ 0, degree(D,x n ) = 0 then 

Proj Xm+lj ... tXrl Zero(AS/JD) = Proj Xm+u ... , Xri _ 1 Zero(AS' / J’ (I S D)) 

3. if degree(A s ,x n ) = 0, degree(D,x n ) / 0 Let D = J2i=o D i x hi A e 
K[x i, • • • ,x n -i], then 

Proj Xm+u ... iXn Z ero(AS/ J D) = |J Proj Xm+u ... iXn _ 1 Zero(AS' / J'(/ S A)) 

i 

4. If degree(A s , x n ) = 0, degree(D,x n ) = 0, then 



Proj Xm+u ... , Xrl Zero{AS/ JD) = Proj Xm+lj ... tXn _ 1 Zero(AS/ JD)) 
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Proof. (1)(2)(4) are obviously according to the definition of projection and the 
fundamental theorem of algebra. For (3), please see [7] for details. 

According to the lemmas given above, we can eliminate one variable. All the 
variables which will be eliminated can be eliminated successively if we compute 
it recursively. 

Let PS be a polynomial set and D be a polynomial in K[ x\, • • • , x n ]. If we 
want to compute Proj Xl: ... jXri Zero(PS/ D). The computing process can be di- 
vided into two steps. The variable ordering should be x n > ■ ■ ■ > x m +i > x m > 

■ ■ ■ > X\. Under this variable ordering, first we will decompose the polynomial set 

PS into a series of ascending sets ASi s.t. Zero(PSfD) = Zero(ASi/ JiD). 

Each ASi is an ascending set for the variable ordering !„>••■> x m +i > x m > 

■ ■ ■ > X\. Then for each Zero(ASi/ JiD) we can compute its projection. In the 
following, we will give the algorithm in detail. 

Step 1. Decompose polynomial set PS into a series of ascending sets ASi s.t. 
Zero{PS/D) = Zero(ASi/ JiD) 

i 

where Ji is the product of the initials of the polynomials in ASi. 

Please see [9] for the detail of the algorithm. 

Step 2. Compute the projection of Zero(AS/ JD), AS is an ascending set. J is 
the product of the initials of the polynomials in AS,D is a polynomial. 

ProjectASCAS, J,D,X,Y) 

Input : 

AS: an ascending set w.r.t X; 

D: a polynomial; 

J: the product of the initials of the polynomials in AS; 

X: a list of all variables with descending order; 

Y: a list of variables to be eliminated, also with descending 
order; 

Output : 

a list, its element is also a~list with the form [’as’j’d 1 ], 
’as 1 is an ascending set and ’ d’ is a'polynomial . 



begin 

y:=First(Y), A:=Last(AS); 

Y’:=Y/y, AS’:=AS/A; 

I : =initial (A) ; 
if Y = [ ] then 

result : = (AS , J*D) ; 

if degree ( A, y) =0 and degree (D,y)=0 then 
result:= ProjectASCAS, J,D,X,Y’) ; 
else if degree (A, y) != 0 and degree (D,Y)=0 then 
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result : = Project AS (AS ’,J’,I*D,X,Y’); 
else if degree (A, y)=0 and degree (D,y) != 0 then 
result : = [] ; 
cf :=coeffs(D,y) ; 
for each c in cf do { 

result:= result union ProjectAS(AS, J,c,X,Y’) ; 

} 

else if degree (A, y) != 0 and degree (D,y) != 0 then 
d:=degree(A,y) ; 

R:=Prem(D~d,A,y) ; 

result : = Project AS (AS 1 , J ’ ,I*R,X,Y) ; 
end if 

return result ; 

end 



In the above algorithm. ”/” means a list gotten by removing the element 
behind it from the list before it; ”!=” means not equal to; initial (A) return 
the initial of A ; degree(A , y) return the degree of the polynomial A w.r.t the 
variable y\ coef fs(D : y) return a list composed of all the coefficients of the poly- 
nomial D w.r.t variable y. Prem(D d , A,y) return the pseudo-remainder of D d 
to A w.r.t variable y. 

4 Application on Geometry Theorem Proving and 
Formula Deduction 

Let K = Q be the rational number field. E = C be the complex number field. 
The corresponding geometric configuration of hypothesis is expressed by a finite 
set of polynomial equations PS = 0 i.e. PS = {pi,--- ,p s },pi = 0, • • • ,p s = 
0,Pi £ Q[x i, • • • , x n ]. The geometric configuration of the conclusion is expressed 
by a polynomial equation C = 0 ,i.e..c = 0,c £ Q[x i,--- ,x n }. We will divide 
the the variables Xi, ■ ■ ■ ,x n into two parts X\, ■ ■ ■ , x m and x m +i • • • , x n . Ap- 
plying Proj Xm+1: ... tXrl to the quasi variety Zero(PS/C), the following theorem 
determines that a geometric theorem is universally true or not. If not, a se- 
ries of polynomial equations and inequations about the variables xi, • ■ ■ , x m are 
given(i.e. Proj Xm+lt ... tXn Zero(PS/C)) , which is the sufficient and necessary 
condition for the theorem to be false. 

Theorem 2. For polynomial set PS and polynomial C as shown above, then 

(1) if Proj Xm+lt ... tXn Zero(PS/C) = 0 and Proj Xm+u ... <Xn Zero(PS/ D) ^ 0, 
then the theorem T is universally true. 

(2) if Proj Xm+lt ... tXn Zero(PS/C) 0, V(ai,--- ,a m ) £ Proj Xm+u ... , Xn Zero(PS 
/C) , then there is (a m +i, • • • , a n ) s.t. PS(ai , • • • , a n ) = 0 and C(a i, • • • , a n ) = 
0; 

(3) if there is (ai, • • • , a n ) s.t. PS(a\, • • • , a n ) = 0 and C(ai, ■ ■ ■ , a n ) ^ 0, then 
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(ai 

5***5 ) G Proj Xm+1 ,... iXn Zero(PS/C), so that Proj Xm+lt ... tXn Zero{PS/C) 

is not empty. 

Proof. (l)we claim that Proj Xm+lr .. <Xn Zero(PS/C) = 0 <f=> Zero(PSfC) = 0. It 
is obviously by the definition of projection. Since Proj Xm+lr .. tXn Zero(PS/C ) = 
0, it follows that Zero(PS/C) = 0, i.e. Zero(PS) C Zero(C) i.e. V a = 
(ai,-- - ,a n ) G E n , pi (a) = 0, • ■ • ,p s (a) = 0 => c(a) = 0 so that the theo- 
rem T = ( PS,C ) is universally true. (2) (3) are obviously by the definition of 
projection. 

Based on the above theorem, we can prove geometric theorem mechanically 
and give the non-degenerate conditions automatically by computing the projec- 
tion of a quasi variety . 

It’s obvious to see that we can predetermine the variables occurring in non- 
degenerate conditions. So it’s convenient for us to observe the range of possible 
value for any variables. As a rule, we will eliminate all the dependent variables . 

Moreover, we can be certain that the conditions found through this method 
are the weakest compared to the condition obtained by Wu’s method and the 
others based on Grobner basis such as Kapur’s and Winkler’s approaches. 

Now we consider the application of projection method on formula deduction. 

Theorem 3. The hypothesis of a geometric statement is expressed by a set of 
polynomial equations PS = 0, then Proj Xm+lt ... tXn Zero(PS) gives a series of 
polynomial equations and inequations which involve the variables x \ , • • ■ , x m 
only. 

While deducing the unknown geometric formula by the projection method, 
all the variables which do not occur in the final formula will be eliminated. 

5 Examples 

In this section, we will give two examples to show that how the projection method 
is applied to automatic theorem proving and formula deduction. The first ex- 
ample is taken from Wang [8]. We will prove it by computing the projection 
of quasi variety and give the non-degenerate conditions and compared the re- 
sult with Wu’s method. The second example is to derive the geometric formula 
automatically from the given hypothesis. 

Example 1. The bisectors of the three angles of an arbitrary triangle, three-to- 
three, intersect at four points. Let the triangle be ZiABC, the two bisectors of 
ZA and ZB intersect at point D. We need to show that CD is the bisector of 
ZC. 

We take the coordinates of the points as A{x\, 0), B(x 2,0), C(0, x^) 1 D{xn, x$). 
The hypothesis of the theorem consists of the following relations. 

H 1 = x$[x\ — (#4 — x\) 2 } — 2x\x^,{xi — x\) = 0 (DA is the bisector of Z CAB) 
H 2 = £3 [x\ — (aq — X 2 ) 2 ] — 2x2X$(x4 — X 2 ) = 0 (DB is the bisector of Z ABC) 
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C 




Fig. 1 . Three bisectors go through the same point 



The conclusion to be proved is C = 0 

C = [xi(x 5 - Xs) + X 3 X4\[X3(X 5 -x 3 ) - 2:22:4] + [x 2 (x 5 ~ *3) +X3X4] [a7 3 (a?5 -x 3 ) - 
*1X4] 

The characteristic set of {Hi,H 2 } is CS = [CijCy 
Ci = x\(x\ — x 2 )x\ + lowerterms 
C 2 = ^3(^1 — x 2 )(xi + x 2 — x<i)x5 + lowerterms 

The pseudo-remainder of the conclusion polynomial C w.r.t the characteris- 
tic set CS is 0. The theorem is true under the non-degenerate conditions which 
are x 3 ^ 0,Xi ^ x 2 , and X\ + x 2 — X4 ^ 0 in Wu’s sense. 

There are three degenerate cases. They are: 

Case 1 ii = x 2 . In this case, A and B are coincide, then ABC is not a real 
triangle anymore. 

Case 2: *3 = 0. In this case, C is on the line AB, the ABC is not a real triangle 
also. 

Case 3:xi + x 2 — x\ = 0. In this case, the intersection point of the bisectors is 
on line x = x\ + x 2 . In Wu’s method, we can’t determine the theorem is true 
or not at this time. If we want to know if the theorem is still true, we should 
put the polynomial P = X\ + x 2 — X4 to the original polynomial set {H\,H 2 } 
to form a new polynomial set {Hi, H 2 , P}, compute its characteristic set again, 
and decide the theorem is true or not under the condition X\ + x 2 — X4 = 0. 

Now we will prove the theorem and give the non-degenerate conditions by 
computing the projection of the quasi variety. 

In this example, £1,2:2, £3 are the free variables and 214,2:5 are the depen- 
dent variables. We compute the projection of Zero({Hi, H 2 } /C) to eliminate 
the variables 2:4,2:51 Proj XitX5 Zero({Hi, H 2 } /C) = Zero({— x 2 + x{\ / {x 3 ,x\ + 
xl})\JZero({-x 2 + xi, x 3 }/{xi}) 
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Fig. 2. The relation among three line segments PQ,PM and PN 



There are only two degenerate cases. They are: 

Case 1: x 2 = x\, x 3 ^ 0. In this case, A is coincidence with B , and C is not on 
the line y = 0. 

Case 2: x\ = X 2 and *3 = 0 and x\ ^ 0, In this case, A is coincide with P,but A 
is not coincide with C . 

This theorem is false only under these two degenerate cases and it is always 
true except these two degenerate cases. 

Example 2. Let R be a point on the circle with diameter AB. At a point P(not A 
or B) of AB a perpendicular is drawn meeting BR at N, AR at M , the circle 
at Q, Find the relation among PQ,PM and PN. 

First, take the coordinate of the points as O = (0,0), A = (xi,0 ),B = 
(x 2 ,0),P = (x 3 ,0),R = (xi,x 5 ),M = {x 3 ,y{),Q = (x 3 ,y 2 ),N = (x 3 ,y 3 ) The 
hypothesis is expressed as following polynomial equations. 



H 1 = x\ + x 2 ( AB is diameter ) 

H 2 = x\ — x% — x\ (|AO| = \RO\) 

H 3 = —y 3 {xi - x 3 ) - (x 5 - y 3 )(x 2 - ^ 3 ) (N, R, B are collinear ) 

= X 5 (x 3 — Xi) — y\(xi — X\) (A, M, R are collinear ) 

H 3 = x\-y\-x\ (|AO| = \QO\) 

Since we want to derive formula about PQ,PM and PN, it is a formula 
about the variables about y 3 ,y 2 ,yi- Since P ^ A, P ^ B and A ^ O, let D = 
{x3 — x\, x 3 — x 2 , Xi} in advance. Then eliminate the variables x 3 , X4, x 3 , x 2 , x\ 
by computing the projection of Zero({Hi, H 2 , H 3 , H 4 , H 5 }/D). 
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Proj X5 ,X4,X3,X2,X\ Zero{{H u H 2 ,H 3 ,Ih,H 5 }/D) 

= Zero({y 3 y i - yQ/{ y 2 }) U Zero({yi}/{y 2 }) U Zero({y 3 } / {y 2 }) 

There are three degenerate cases. They are: 

Case 1: y 3 y\ —y 2 =0,y 2 j^0 i.e. When R doesn’t coincide with A or B , we have 
| PiV | * \PM\ = \PQ\ 2 . 

Case 2: y 3 = 0,y 2 / 0 i.e. R coincides with B 
Case 3: y 3 = 0,y 2 ^ 0 i.e. R coincides with A 

6 Conclusion 

We give an algorithm to compute the projection over an algebraic closed field. 
Applying this algorithm to automatic theorem proving, we can get the weakest 
non-degenerate condition for which the theorem is true. In fact, we can get 
the sufficient and necessary condition for a geometric theorem to be false by 
computing the projection of a quasi variety. This algorithm also can be applied 
to automatic geometric formula deduction. There are more than one hundred 
geometric theorems which have been proved by this method. Experiments show 
that the projection of the quasi variety can be computed out if we can get the 
zero decomposition for the given polynomial system. 
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Abstract. In this paper we present a classification of 3-revolute-jointed 
manipulators based on the cuspidal behaviour. It was shown in a pre- 
vious work [16] that this ability to change posture without meeting a 
singularity is equivalent to the existence of a point in the workspace, 
such that a polynomial of degree four depending on the parameters of 
the manipulator and on the cartesian coordinates of the effector has a 
triple root. 

More precisely, from a partition of the parameters’space, such that in 
any connected component of this partition the number of triple roots 
is constant, we need to compute one sample point by cell, in order to 
have a full description, in terms of cuspidality, of the different possible 
configurations. 

This kind of work can be divided into two parts. First of all, thanks 
to Groebner Bases computations, the goal is to obtain an algebraic set 
in the parameters’space describing the cuspidality behavior and then to 
compute a CAD adapted to this set. 

In order to simplify the problem, we use strongly the fact that a manip- 
ulator cannot be constructed with exact parameters, in other words, we 
are just interested in the generic solutions of our problem. 

This consideration leads us to work with triangular sets rather than with 
the global Groebner Bases and to adapt the CAD of Collins as we will 
just take care of the cells of maximal dimension. 

Keywords: Groebner Basis, Cylindrical Algebraic Decomposition, Tri- 
angular Sets, Serial Manipulators, Semialgebraic Sets 



1 Introduction 

Industrial robotic 3-D OF manipulators are currently designed with very simple 
geometric rules on the designed parameters, the ratios between them are always 
of the same kind. In order to enlarge the possibilities of such manipulators, it 
may be interesting to relax the constraints on the parameters. 

The behavior of the manipulators when changing posture depends strongly on 
the design parameters and it can be very different from the one of manipulators 
commonly used in Industry. 



F. Winkler (Ed.): ADG 2002, LNAI 2930, pp. 31-43, 2004. 
(c) Springer- Verlag Berlin Heidelberg 2004 
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P. Wenger and J. El Omri [9], [16] have shown that for some choices of the 
parameters, 3-DOF manipulators may be able to change posture without meeting 
a singularity in the joint space. This kind of manipulators is called cuspidal. 
It is worth noting that in case of obstructed environment, this property would 
yield more flexibility which can be very useful in practice for industrial purpose. 
They succeed in characterizing 3-revolute jointed manipulators using a homotopy 
based classification scheme [15], but they needed general conditions on the design 
parameters, more precisely they wanted to find answers to the following issues: 

— Problem 1 : For given parameters, is the manipulator cuspidal? 

— Problem 2 : How to obtain an algebraic set in the parameters ’space such 
that in each connected components of the complementary, the cuspidality 
behavior is fixed? And then can we caracterize each connected component 
with a test manipulator? 

We restrict the study to 3-DOF manipulators as described in figure 1. The 
geometrical modelisation of the manipulator is explained in the first section of 
this paper. As recalled in the second section, testing if such a manipulator is 
cuspidal or not is equivalent to deciding if an algebraic set has real roots or not, 
this set corresponding to the existence of a triple root of a given polynomial of 
the fourth degree. 

We manage to answer both questions under a few hypotheses which are not 
restricting ones according to roboticians : for example it is impossible to con- 
struct, in practice, a manipulator whose parameters lies on a strict hypersurface 
of the space of the parameters. 

We explain in the third section how we obtained the algebraic set, depending 
only on the parameters, which complementary interests us for the study of the 
cuspidality behavior. Then in the fourth section, we show how to compute a par- 
tition of the space of parameters, such that in each cell of maximal dimension, 
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the behavior of the manipulator is known (cuspidal or not), while the other cells 
are embedded inside a strict algebraic subset of the parameters’space. 

Many authors propose to use computer algebra tools for solving geometrical 
models of serial manipulators (see [7] for example), mainly using resultants the- 
ory. In this paper, the main computational objects are Grobner bases (as in [4]) 
and triangular sets ([2], [14]). 

2 The Kinematic Map 

The 4 design parameters are g^, ds, and r^. In order to normalize 
(g? 2 , d 3 , g?4, T2) and to reduce the number of parameters, we assume that d'2 is 
equal to 1. 

So the space of design parameters is (R + ) 3 . Along our study of the behavior 
of the manipulator, we will work into two spaces: 

— the joint space described by the joint variables (0i, 02 , 0 3 ), this space is 
isomorphic to ] — 7r,7r] J as the joint limits will be ignored. Moreover, as the 
manipulator turns around the axe of the first revolute joint, the joint space 
can be represented by the torus (in the variables (02, 0 3 ))- 

— the task space representing the position of the end-effector in the Cartesian 
coordinates (x,y,z), which is isomorphic to K 3 if we suppose there are no 
obstacles. 

The kinematic map / maps the joint space on the task space : 

/ : ] -7T,7r] 3 > M 3 

(01,02,03) 1 — » (x,y,z) 

The image of this map in the task space is called the workspace. 

To express /, we just have to express the change of bases between the bass of 
the origin and the basis of the effector, which gives us the following expression 
for f: 



! x = (g? 3 + cos 03(4) (cos 0i cos 02) + (r2 + ch sin 0 3 ) sin 0i + cos 0i 
V = (^3 + cos 0 3 d4)(sin 0i cos 02) — (7-2 + d± sin 0 3 ) cos 0i + sin 0i (1) 

z = (d 3 + cos 0 3 di) sin 0 2 

For the user, it is of first importance to know the workspace, he wants to 
be able to divide the task space into accessible or not accessible zones, and to 
know the number of joint’s configurations corresponding to the same position. 



3 Algebraic Characterization of Cuspidal Manipulators 

In the classical cases used in Industry, manipulators, to change posture, need 
to pass through a singularity of the joint space, in other words the end-effector 
must bump into the frontiere of the workspace. 
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But this behavior is not general at all. It was shown by P. Wenger that 
a 3-DOF manipulator can execute a non singular change of posture if and only 
if there exist at least one point in its workspace with exactly three coincident 
inverse kinematic solutions (corresponding in a cross section of the workspace to 
cusp points, hence the word cuspidal). 

As it is difficult to express such a condition directly from the kinematic map, 
the idea is to eliminate joint variables from the system, let say 6 \ and 6 2 , in 
order to obain a condition over the last joint variable. 

By taking Si = sin(9i ) and Ci = cos(9i ), i = 1 ... 3, and adding the algebraic 
relations s 2 + cf = 1 i = 1 . . . 3 to the equations defining /, we have to study an 
algebraic system of equations. 

In order to eliminate 9\ and 62 , we compute a Groebner Basis of the system 
under the following elimination order: 



[ci,si,c 2 ,s 2 ] > [C3, s 3 , r 2 , d 3 , d 4 , y, z]. 

The Groebner basis G 3 of the Elimination Ideal is composed of two polyno- 
mials : 

cl + sl -1 

and 

TO5C3 + TO4S3 + TO3C3S3 + m 2 c 3 + 7711S3 + m 0 

where : 

ry T) 1 2 1 ( -F2— |— 1 — L 

m 0 = Z - R + r 2 -\- 

mi = 2r2a?4 + (L — R — l)d4r 2 

m 2 = (L — R — l)d 4 d 3 

m 3 = 2 r 2 d 3 dl 

m 4 = d\(rl + 1) 

. ms = d\d\ 

with R = x 2 + y 2 + z 2 , Z = z 2 and L = d\ + d 2 + r 2 . 

After a change of variables in t = tan(^-), the second polynomial becomes: 

P(t) = at 4 + bt 3 + ct 2 + dt + e 



with: 

( o = ms - m2 + mo 
b = —2 m 3 + 2mi 
c = —2ms + 4rri4 + 2mo 
d = 2 m 3 + 2mi 
_ e = ms + m2 + mo 

Deciding if a manipulator is cuspidal is equivalent to deciding if this poly- 
nomial P of degree 4 in t (whose coefficients are polynomial with respect to 
x, y , z, g? 3 , r 2 ) admits real triple roots. 

It is important to note that every solution in t = tan(^) is uniquely lifted in 
a 3-uplet (0i, 62 , $ 3 ), except ii z = 0 or x 2 + y 2 = 0, but this case will be treated 
later. We will show that the design parameters 03 ,^ 4 , r 2 of manipulators such 
that z(x 2 + y 2 ) = 0 are in a strict hypersurface of R 3 . So, finding parameters’ 
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values defining cuspidal manipulators remains to find the values of cLt, cfe, r 2 such 
that P(t) has triple points by solving the following system (as the boundaries 
of the workspace form a revolution surface around the axis Oz, it is generically 
zero-dimensional once the parameters are fixed) : 



P =0 

dP 
Qt 
dPP 
dt 2 



= 0 
= 0 



( 2 ) 



4 The Algebraic Set Defining the Classification 

We want to obtain a partition of the space of parameters, such that the cuspi- 
dality behavior of the manipulators corresponding to the same cell is identical. 

As it is impossible to construct in practice a manipulator whose parameters 
correspond exactly to given parameters, we just deal with “generic” solutions. 
In other words, we are only interested in the cells of maximal dimension. In each 
of those cells we want that the number of real roots of our system is constant 
(in the generic cases it appears that the dimension of thr system when we fix 
the parameters is zero). 

In order to obtain some conditions on the parameters, we need to eliminate 
the variables t, Z 1 R from the system of equations by computing an eliminating 
polynomial depending on the parameters d 4 ,d 3 ,r 2 , and on one of the three 
variables t, R and Z. But, moreover, we need to caracterize the solutions de- 
pending on these conditions which impose to compute simple expressions of the 
coordinates of the solutions with respect to the two remaining variables. 

The eliminating polynomial could be computed using Grobner basis as well 
as resultants, but none of the two algorithms may give simple expressions of the 
solutions as expected. 

For example, we have computed a Grobner Basis of the system with respect 
to an elimination order. This is possible using the most recents algorithms from 
J.C.Faugere [5]. But the result obtained is very huge and difficult to use because 
all the multiplicities as well as non real components of the variety, are kept along 
the computation. 

So, we choose to represent the solutions of system (2) as regular zeroes of 
regular and separable triangulars sets (with respect to the terminology of [2]). 
This ensure that the non “generic” solutions correspond to values of the pa- 
rameters contained in strict hypersurfaces of the space of the parameters, that 
we will work with equidimensionnal ideals and that “generic” specializations of 
the parameters will lead to zero-dimensional systems without multiple complex 
roots (critical for the study of real roots) . 

As the direct computation do not give a pleasant result, we decided to modify 
the variables and to simplify the problem. Thanks to the following change of 
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variables: 



u = 


e — a 


— 2 (d\ + d\ + r \ — 


Xi = d — 


CO 

1 

T3 

II 

-O 


hU = Sr 2 d 3 d\ 


X 2 = 


U + 2 a — c 


= Adl(dl-rl-l) 


h = 


T2/d 3 




x 3 = 


b-hU 


= Ar 2 di{l — d 3 d±) 



R - l)d^r 2 



the Grobner basis of the new system with respect to the lexicographic order : 



t > a > U > X! > X 2 > -X 3 > h 

is so easy to compute that it can be done very quickly with the algorithm of 
Buchberger implemented under MAPLE. 

The basis obtain as the following form : 

' surf u (U,X 1 ,X 2 ,X 3 ,h) =0 
5 a,i(o, U, Xi, X 2 , X 3 , h) = 0 
9 a, 2(0, U, X\, X 2 , X 3 , h) = 0 



< ffa,14(o, u, X±, X 2 , A3, h) = 0 
a, U , Xi,X 2 , X 3 , h) = 0 
a, U , Xi,X 2 ,X 3 , h) = 0 



„ 9t,ia(t , a j y j Ai,A 2 , X 3, h) — 0 

where : 

— the polynomial is of degree 1 in the variable a. 

— the polynomial g t ,i is of degree 1 in the variable t. 

We do not give here the full polynomials because of their size (more than 320 
monomials for each of them) . We extract from the system the following triangular 
SGt* 

f surf u {U,X 1 ,X 2 ,X 3 ,h) = 0 

{ g a ,i = lc a (X 1 ,X 2 ,X 3 ,h)a + c a (U,X 1 ,X 2 ,X 3 ,h) = 0 (3) 

[ 9t,i = lct(X i,X 2 , X 3 , h)t + ct(U,a , Xi, X 2: X 3 , h) = 0 

This triangular set T extracted from the system is regular and separable, it 
means that the product of the leading coefficients of T (lctlc a lc sur f v ) as well 
as the derivatives of g a , 1 and g tt 1 with respect to their main variables, do not 
vanish on any prime component of the ideal < T > (this has been verified by 
localizing the system). 

Hence the saturation ideal of < T > is an equidimensional ideal which ensures 
us that the different projections (eliminations of variables) we will apply on this 
ideal will have a good behaviour (good dimension of the varieties obtained at 
each step). In fact we are working with the Zariski closure of the variety of 
the regular zeros of < T > which is here an equidimensional envelope of the 
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variety of the zeros of the components of dimension 3 of the system (2). Note 
that basically, the system is of dimension 4, but one can easily show that the 
components of dimension 4 are included in the hypersurface = 0. 

So coming back to our first variables, we can define the generic solutions of 
the problem as regular roots of a triangular system with the following shape: 

{ surf(R,d 4 ,d 3 ,r 2 ) = 0 

lc z (d 4 ,d 3 ,r 2 )Z + tr z (R,d 4 ,d 3 ,r 2 ) = 0 (4) 

lct(d 4 , d 3 ,r 2 )t + tr t {R, Z 1 d 4 ,d 3 ,r 2 ) = 0 

Now to obtain a description of the complementary of our algebraic set, we 
can compute directly a Cylindrical Algebraic Decomposition adapted to this set 
and then just keep the cells of maximal dimension. But the computation of the 
discriminant of the polynomial surf would be a limitating step. 

Let’s study the system more deeply. 

The set Vi c (Z,t) = {(d 3 ,d 4 ,r 2 )\lc z (d 3 ,d 4 ,r 2 )lc t (d 3 ,d 4 ,r 2 ) = 0} defines an 
hypersurface in the parameters’space, which is a strict closed set. We can remark 
that lc z = 0 has no real roots, so VJ C = {(cfe, c? 4 , r 2 ) \ lct(d 3 , d 4 , r 2 ) = 0}. Studying 
the generic solutions of the first system is exactly avoiding the parameters on 
this Vic and then parameters for which the number of real roots of surf(R) 
varies. It means the parameters lying in the hypersurface : 

{(d 3 ,d 4 ,r 2 )\lc sur f = 0 or discrim(sur f ', R) = 0} 

The real hypersurface lc sur f = 0 is empty and the discriminant factorise in 
three factors: 

- fi = d 3 - d\ + r| 

— f 2 = d 2 d® — + 3d\d\r 2 — 2d\d\ + 2d\d 2 — 2r 2 d 2 df + d\d 2 + 3 d\d\r% — 

d^r 2 — 2 r 2 d\ — — d\ + r%d\ + d 2 r% + 2 r%dl 

— f 3 — 2 r 2 d 2 d 3 + 2 r%d 4 d 3 + 6 r 2 d 4 d 3 + 6 r 2 d 2 d 3 + ‘ir 2 0 d 4 + 6 r 2 0 d 3 d 4 + 15r2d 3 d4 + 
ISrfdgd! + 12 r 2 <i§d 4 — Ar 2 d\d\ + 6 r 2 d 4 — Ad 3 d\ + 6 d 3 df + r\ 2 d\ — 4r 2 d 2 d 3 + 
Ar^d^d 2 + 20 r 2 d 3 <i| — 8 r 2 d\d\ + 8 r 2 d 3 d\ + d\r 2 + 2 d 2 d 3 r 2 + 2 d\d 3 r\ + d\d\ — 
1 2r| d\ — 2r| d\ + d^r^ + 4r® d 2 d\ — 4r \ d 2 d\ + 4r 2 e?| d\ + 4r| c?| d\ + 4r ® d\ + 
dffd\ — 4d 3 °d4 + Qr^df^d'l 

The real hypersurface f 3 = 0 is also empty from the real point of view. Also 
we must remember that the variables Z and R— Z are strictly positive. Substi- 
tuting Z = 0 or R — i? = 0 in the first system permits us to find two polynomial 
conditions on parameters, let say Z 3 (d 3 , d 4 , r 2 ) = 0 and RZ 3 (d 3 ,d 4 ,r 2 ) = 0. 

The polynomial RZq, does not present any real root. The polynomial Zq 
factorises in two polynomials : 

— 2T 0 i = d 3 r% + d 2 — 2 d 3 + d 3 — d 2 + 2 d 3 d 2 — d 2 d 2 
— Z 32 = d 2 r 2 + d 2 + 2d\ + d\ — d\ — 2d 3 d\ — d 2 d 2 

The algebraic set to avoid : 



{Zqi — 0, Z 32 — 0, /i — 0, f 2 — 0, let — 0} 
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gives us a partition of the space of parameters (where d 3 >0,^4 > 0,r2 > 0), 
such that in each cell of maximal dimension the number of cusp points in a cross 
section of the workspace is constant. 

5 Partition’s Cells Computation 

A way to represent such cells is now to compute a Cylindrical Algebraic Decom- 
position (CAD - see [ 3 ]) of R 3 adapted to this set of polynomials. But it may 
give us a result very huge and difficult to analyse in practice, and with lots of 
cells we are not interested in. 

In order to obtain a decomposition easier to manipulate, we use the fact 
that we are just interested in the “generic” solutions, so in the cells of maximal 
dimension in our partition of the space of parameters. 

So we adapt the Collins ’Algorithm to our case. 

— Projection phase 

At each step we keep the leading coefficients of the polynomials, the resul- 
tants and the discriminants, keeping in mind that the polynomials which 
have no real roots do not interest us. As we are not interested in the mul- 
tiplicities of the singularities and intersections between hypersurfaces, we 
won’t compute the subresultants. 

— Lifting phase 

As we are just interested in the cells of maximal dimension, we do not need 
to compute real algebraic numbers and we can work only with rationnal test 
points (so we won’t have to work with towers of extention for example, or 
other ways of coding real algebraic numbers). 

5.1 The Projection Phase 

We want to give a description of the connected components of maximal dimen- 
sion of the complementary of the algebraic set A = {/1, /2, , let , Z01, ^02} where: 

' fi = d\- d\ + rl 

/2 = d|d 3 — d\d\ + Zdjdfr?; — 2e?|d| + 2 d\d\ — 2 r\d\d\ + d\d\ + Hd\d\r^ 

— d\r\ — 2 r%df — — d\ + + d\rl + 2 r^d\ 

Z m = c?|r| + — 2 dl + df — d\ + 2 d$d\ — d\d\ 

Z02 = + 2 — d\ — 2 d^d\ — d^d^ 

k lc t = - d 3 + d^l + C?4 

Let 

Pl?OJ d4 (A) = {AnQ[d 3 ,r 2 ],d*sc(P, d 4 ),res(P, Q,d 4 ),k(P) | P,Q<=A}. 

It is the smallest family of polynomials in Q[d 3 , r 3 ] which keep in mind the struc- 
ture of the connected components of maximal dimension of the complementary 
of A. After computations we obtain: 



PROJ di (A) = {C 1 ,...,C 6 }. 
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where: 

d 3 — 1, 

— 1 + 2 d 3 , 

-da 2 + d 3 4 + 2 d 3 2 r 2 2 + r 2 2 + r 2 4 , 
d 3 r 2 +2 d 3 — 2 d 3 r 2 — 4 d 3 -f- 3 r 2 T 3 T r 2 , 

-3 d 3 2 + 3 r 2 2 + d 3 4 r 2 2 + 2 d 3 4 + 3 d 3 2 r 2 4 + 4 d 3 2 r 2 2 + 3 r 2 6 + 6 r 2 4 , 

3 d 3 2 - 3 r 2 2 + 10 d 3 4 - 10 d 3 2 r 2 2 - 3 r 2 4 + 10 d 3 r 2 2 - 10 d 3 3 - 5 d 3 6 - 
11 d 3 4 r 2 2 — 7 d 3 2 r 2 4 — r 2 6 + 8 d 3 3 r 2 2 + 8 d 3 r 2 4 + 2 d 3 r 2 6 + 6 d 3 3 r 2 4 + 

6 d 3 s r 2 2 + 2 d 3 7 . 

We can visualise this algebraic set and the partition of the parameters ’space 
(d, 3 ,r 2 ) it implies, on the following figure. 

Then we define by the same way 

PROJ r2 {PROJ di (A)) = {-1 + 2 d 3 ,d 3 - 1,2 d 3 - 3, -3 + 2d 2 }. 



'Ci = 
C 2 = 
C 3 = 
C 4 = 
C 5 = 

c 6 = 



5.2 The Lifting Phase 

At each step of the lifting phase we use the variant of Descarte’s rule of sign based 
method proposed in [12] (many other methods like [6], [8] or [1] should be used) 
to find isolating intervals of the real roots appearing along the computation. 

For the first step, given the list of polynomials in the parameter d 3 , we isolate 
in intervals with rational endpoints all the real roots of these polynomials. Then 
we choose the test values for d 3 to be the middle points of each interval between 
the isolating intervals. We obtained the following test values for d 3 : 
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- 4 1} = 

7(2) _ 



219 

1024 



e]0,Pi 



= - el PuP 2 



4 3) = 



4 
1139 



1024 
(4) _ 2791 



4 5) = 



2048 

1011 

512 



elP 2 ,P 3 



e]i^P 4 [, 

e]P 4 ,oo[, 



After substituting those values in the list {Ci, Ce}, we obtain 6 polynomials 
in r 2 . In each case we isolate the real roots and obtain the following test values: 



for 4^ = ^ , 

3 1024 



219 , r 403 

: the test values for r 2 are { 



403 827 215 13657 



for 4 2) = - 
3 4 



{ 



339 371 939 67189 



— for 4 3 ) = 



for = 



2048’ 1024’ 2048’ 128000 
1139 r 259 785 32907, 

1024 : 3 2048’ 2048’ 64000 
2791 r 181 45439 , 

l -1 nr, a J 1 oonnn J , 



4096’ 2048’ 4096’ 1024’ 64000 

}, 






and for d y 3 = 



2048 ’ L 1024’ 128000 

(5) 1011 



512 



{!}• 



Then we do the same in the set of polynomials A. Over each 15 cells of the 
plane (r 2 ,d 3 ) (discribed in figure 2) we obtain 7 test values in the parameter d 4 
(because the algebraic set A is generically composed of 6 leaves as the hypersur- 
face / 2 = 0 is composed of two leaves) , so we eventually found 105 test points 
(each one corresponding to a cell of maximal dimension). 



5.3 Results 

For each sample point, after specialisation in the given parameters, it is then 
sufficent to solve the following zero-dimensionnal system (counting the number 
of real roots) to get the number of cusp points in a cross section of the workspace 
corresponding to the selected parameters: 

P =0 

f =° 

< 1^=0 ( 6 ) 

T 3 Z - 1 =0 

, T*(R -Z)- 1 = 0 

where the equations T 3 Z — 1 = 0 and T^{R — Z) — 1 = 0 discriminate the 

admissible solutions. 

This can be done computing a Rationnal Univariate Representation of the 
system usign algorithms described in [10], [11] (other strategies like [13] should 
also be applied without any trouble to compute such a rational parametrization) 
and then isolating the real roots of this RUR with the algorithm described in [12] 
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(this last computation is not the main one : numerous variants [6], [8] or [1] could 
be used as well as Sturm sequences).; 

As Z = z 2 , the number of cusp points appearing in a cross section of the 
workspace, equals twice the number of real roots of the RUR. 

Here are the results for each sample point: 



(d 3 ,r 2 ) \ c?4 


1 


2 


3 


4 


5 


6 
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(1,1) 


0 


0 


4 


4 


2 


0 


0 


(1,2) 


0 


4 


4 


4 


2 


0 


0 


(1,3) 


0 


4 


4 


4 


2 


0 


0 


(1,4) 


0 


4 
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(1,5) 
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0 
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(2,1) 
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(2,2) 
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0 


(2,3) 


0 


4 


4 


4 


2 


2 


0 


(2,4) 


0 


4 


4 


2 


2 


2 
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(3,1) 


0 


4 


4 


4 


2 


2 


4 


(3,2) 


0 


4 


4 


4 


2 


2 


4 


(3,3) 


0 


4 


4 


2 


2 


2 


4 


(4,1) 


0 


4 


4 


4 


2 


2 


4 


(4,2) 


0 


4 


4 


2 


2 


2 


4 


(5,1) 


0 


4 


4 


2 


2 


2 


4 



In each line (*, j) of this table you can read the number of cusp points appear- 
ing in a cross section of the workspaces of the seven test robots above the cell 
(i,j) (see figure 2) corresponding to the seven distinct test values of obtained 
in the fiber over the point (dg^r^)- 



Example 

Here are the different cross sections of the workspace obtained for the sample 
points from the fiber above the cell (5, 1). The drawings are from SCILAB. The 
value of the parameter c ?4 goes increasing from left to right. 

We can verify on those drawings that the number of cusp points corresponds 
to the corresponding number in the table. 



(d 3 ,r 2 ) \ e?4 


1 2 


541 


6 7 


(5,1) 


0 4 z 


12 1 


>24 



Now, if, given parameters, the user wants to know in which cell lies the ma- 
nipulator, he can do it using the lifting phase. 

The first step of the search is to locate the parameter d 3 between the real 
roots of the set PROJ r2 (PROJd 4 (A)) = {—1 + 2d 3 , d 3 — 1, 2 d 3 — 3, — 3 + 2 d§}. 
As each root is given by an isolating interval with rationnal bounds, which can 
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( 5 , 1 , 1 ) ( 5 , 1 , 2 ) ( 5 , 1 , 3 ) 



0 



Fig. 3. Fiber 

be refined if necessary, it is easy to find the first indice of the cell (the one 
corresponding to g^). 

Then we solve the differents equations of the set PROJd 4 (A) = {Ci, Cq}. 
evaluated in the given parameter (so we just have to find the real roots of 
low-degree polynomials in one variable). By the same way than before we locate 
the parameter r 2 and obtain the second indice. By repeating this process with 
these two parameters and the set A we obtain the triple of indices corresponding 
to the cell in which lies the given robot. 

Before concluding this paper, we must remark that the classification we have 
made is not optimal as some neighbouring cells correspond to the same number 
of cusp points. We may explain that by the fact that some leaves of our hypersur- 
faces in A, do not give us informations on the real roots we are searching for. We 
are now trying to find a minimal decomposition in cells of maximal dimension 
of the space of parameters. 

Conclusion 

In this paper, we manage to compute an algebraic set depending on the design 
parameters, such that in each connected component of its complementary the 
cuspidality behavior is the same . Thanks to an adapted CAD, we obtained 
a decomposition of this set in cells and at least one sample point, with rationnal 
coordinates, in the interior of each cell of higher dimension. 

Among those sample points, P. Wenger’s team found just one class of ho- 
motopy, in the torus representing the joint space, over the seven classes they 
encountered before for the cuspidal manipulators [15]. 

Generalization of this work to manipulators with one more design parameter 
(the parameter of the DH parameters) would permit to find representants of 
other classes. We are now working on. 

Future research work is to try to generalize our study to other problems 
involving system of polynomial equations depending on parameters. 
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Abstract. We introduce a software package, MMP/Geometer, devel- 
oped by us to automate some of the basic geometric activities including 
geometric theorem proving, geometric theorem discovering, and geomet- 
ric diagram generation. As a theorem prover, MMP/Geometer imple- 
ments Wu’s method for Euclidean and differential geometries, the area 
method and the geometric deductive database method. With these meth- 
ods, we can not only prove difficult geometric theorems but also discover 
new theorems and generate short and readable proofs. As a geomet- 
ric diagram editor, MMP/Geometer is an intelligent dynamic geometric 
software tool which may be used to input and manipulate geometric di- 
agrams conveniently and interactively by combining the idea of dynamic 
geometry and methods of automated diagram generation. 

Keywords: Geometry software, automated reasoning, geometric the- 
orem proving, geometric theorem discovering, geometric diagram gener- 
ation, intelligent dynamic geometry. 



1 Introduction 

MMP (Mathematics-Mechanization Platform) is a stand-alone software platform 
under development [12], which is supported by a Chinese National Key Ba- 
sic Research Project “Mathematics Mechanization and Platform for Automated 
Reasoning.” The aim of MMP is to mechanize some of the basic mathematical 
activities including automated solution of algebraic and differential equations 
and automated geometric theorem proving and discovering, following Wu’s idea 
of mathematics mechanization [33, 36]. MMP is based on algorithms for symbolic 
computation and automated reasoning, and in particular the Wu-Ritt charac- 
teristic set (CS) method [36, 30]. MMP also implements application packages 
in automated geometric reasoning, automated geometric diagram generation, 
differential equation solving, robotics, mechanism design and CAGD. 

* This work was supported by a National Key Basic Research Project (NO. 
G19980306) and by a USA NSF grant CCR-0201253. 



F. Winkler (Ed.): ADG 2002, LNAI 2930, pp. 44-66, 2004. 
(c) Springer- Verlag Berlin Heidelberg 2004 
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MMP/Geometer is a package of MMP for automated geometric reasoning. 
The aim of MMP / Geometer is try to automate some of the basic geometric ac- 
tivities including geometric theorem proving, geometric theorem discovering, and 
geometric diagram generation. The current version is mainly for plane Euclidean 
geometries and the differential geometry of space curves. 

The goal of MMP / Geometer is to provide a convenient and powerful tool to 
learn and use geometry by combining the methods of geometric theorem proving 
and geometric diagram generation. The introduction of computer into geometry 
may give new life into the learning and study of the classical field [10]. Geometry 
is at the heart for many engineering problems from robotics, CAD, and computer 
vision. We expect that MMP/Geometer may have applications in these fields. 
Actually, some of the methods implemented in MMP/Geometer are directly 
targeted at engineering problems [11, 15, 16, 17]. 

1.1 Automated Geometric Theorem Proving and Discovering 

Study of automated geometric theorem proving (AGTP) may be traced back to 
the landmark work by Gelernter and his collaborators [18] in the late 1950s. The 
extensive study of AGTP in the past twenty years is due to the introduction of 
Wu’s method in late 1970s [33], which is surprisingly efficient for proving difficult 
geometric theorems. Inspired by this work, many successful methods of AGTP 
were invented in the past twenty years. For a recent survey of these methods, 
please consult [7]. 

AGTP is one of the successful fields of automated reasoning. There are 
few areas for which one can claim that machine proofs are superior to human 
proofs. Geometry theorem proving is such an area. Our experiments show that 
MMP/Geometer is quite efficient in proving geometric theorems. Within its do- 
main, it invites comparison with the best of human geometric provers. Precisely 
speaking, we have implemented the following methods. 

Wu’s method might be the most powerful method in terms of proving difficult 
geometric theorems and applying to more geometries [33, 36]. Wu’s method 
is a coordinate-based method. It first transfers geometric conditions into 
polynomial or differential equations in the coordinates of the points involved, 
then deals with the equations with the characteristic set method. 

The area method uses high-level geometric lemmas about geometric invari- 
ants such as the area and the Pythagorean difference as the basic tool of 
proving geometric theorems [8] . The method can be used to produce human- 
readable proofs for geometric theorems. 

The deductive database method is based on the theory of deductive 
database. We may use it to generate the fixpoint for a given geometric config- 
uration under a fixed set of geometric rules or axioms [9]. With this method, 
we can not only find a large portion of the well-known facts about a given 
configuration, but also to produce proofs in traditional style. 

For almost every method of AGTP, there is a prover. We will not give de- 
tailed introduction to existing geometric software packages. A survey may be 
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found in [21]. Comparing with previous provers, MMP/Geometer has the follow- 
ing distinct features. First, it implements some of the representative methods for 
AGTP, while most previous provers are for one method. An exception is Geom- 
etry Expert (GEX), which also implements the methods mentioned above [13]. 
The difference is that GEX only implements a simple version of Wu’s method 
due to the lack of an implementation of Wu-Ritt’s zero decomposition theo- 
rem. Also, we implements general methods of automated diagram generation in 
MMP/Geometer and GEX can handle diagrams that can be drawn with ruler 
and compass constructions. Another reason for us to develop a new software 
package is that GEX is developed under Borland C++ which is not available 
anymore. Second, MMP/Geometer is stand-alone, while most of the previous 
provers are implemented in Lisp or Maple. Third, MMP/Geometer is capa- 
ble of producing human-readable proofs and proofs in traditional style. Finally, 
MMP/Geometer has a powerful graphic interface, which will be introduced in 
the next subsection. 

1.2 Automated Geometric Diagram Generation (AGDG) 

It is often said that a picture is more than one thousand words. But in reality, it is 
still very difficult to generate large scale geometric diagrams like those in Fig. 10 
with a computer. With MMP/Geometer, we intend to provide an intelligent tool 
for generating and manipulating such diagrams. By implementing general AGDG 
methods [11, 24, 20], MMP/Geometer may be used as a general diagram editor. 
Also, AGDG methods have direct applications in parametric CAD [20, 17, 15], 
linkage design [16], etc. 

For AGTP, the diagram editor may provide a nice graphic user interface 
(GUI). Also, AGDG methods may enhance the proving scope for AGTP methods 
with constructive statements as input by finding the construction sequence for 
problems whose ruler and compass construction procedure is difficult to find, as 
shown by the example in Fig. 9. 

Dynamic geometry software systems, noticeably, Gabri [23], Geometer’s 
Sketchpad [22], Cinderella [28], and Geometry Expert [13] may generate dia- 
grams interactively based on ruler and compass construction. These systems are 
mainly used to education and simulation of linkages. As compared with models 
built with real materials, visual models built with dynamic geometry software 
are more flexible, powerful, and more open for manipulation. 

It is well known that the drawing scope of ruler and compass construction 
has limitations. To draw more complicated diagrams, we need the method of 
automated geometric diagram generation, which has been studied in the CAD 
community under a different name: geometric constraint solving (GCS) and with 
a different perspective: engineering diagram drawing. GCS is the central topic in 
much of the current work of developing intelligent and parametric CAD systems 
and interactive constraint-based graphic systems [11, 24, 20, 19]. 

In MMP/Geometer, by combining the idea of dynamic geometry and AGDG 
we obtain what we called the intelligent dynamic geometry , which can be used to 
input and manipulate diagrams more conveniently. It can be used to manipulate 
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Fig. 1 . Simson’s Theorem 



geometric diagrams interactively as dynamic geometry software and does not 
have the limitation of ruler and compass construction. 

The rest of the paper is organized as follows. Section 2 describes the input for- 
mats of a geometric statement. Section 3 reports the implementation of methods 
of AGTP. Section 4 reports the implementation of the methods of AGDG. 

2 Input of a Geometric Statement 

The user may input a geometric statement to MMP/Geometer in two ways: 
either by typing the English description of the statement or by drawing the 
diagram of the statement on a computer screen with a mouse. After a statement 
is inputted, it may be described in four forms: the natural language form , the 
constructive form , the predicate form , and the algebraic form. The purpose of 
using different input forms is that each input form has its merit. For instance, the 
natural language input may be the favorable choice for high school students. The 
constructive form is the input form for several proving methods, like the area 
method [8]. The predicate form is the most general way of describing a statement. 
We use Simson’s Theorem to illustrate these forms. 

Natural Language. MMP/Geometer accepts geometric statements described 
with a precisely defined pseudo-natural language. Detailed definition of this 
language may be found in [12]. MMP/Geometer may convert a statement in 
natural language into all other forms. Using this language, Simson’s Theorem 
may be described as follows. 

geom(" Example Simson. Let D be a~point on the circumcircle 
0 of triangle ABC. E is the foot from point D to line AB . 

F is the foot from point D to line BC. G is the foot from 
point D to line AC. Show that points E, F, and G are 
collinear . " ) ; 

In the above example, geom is an MMP/Geometer command which accepts 
a geometric statement in any form and converts it to all other possible forms. 
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Constructive Form. In this form, the geometric objects in the statements can 
be described by a sequence of geometric constructions. A construction is to 
generate a new geometric object with lines and circles. Detailed description 
of the constructions used in MMP/Geometer may be found in [2, 5, 12]. 
A statement in constructive form may be converted to all other forms. The 
constructive form of Simson’s Theorem is as follows. 

geom( [ [ [P0INT,A,B,C] , [CIRCUMCENTER,0,A,B,C] , [0N,D, [CIR,0,A]] , 
[FOOT , E ,D , A , B] , [FOOT ,F , D , B , C] , [FOOT , G , D , A , C] ] , [ [coll ,E,F , G] ] ] 

); 

Generally speaking, a statement in constructive form is a pair [cs, c], where cs 
is a construction sequence used to generate the geometric objects in a state- 
ment and c is the set of conclusions. 

Predicate Form. This is a natural way to describe a geometric statement. The 
hypotheses and conclusions are represented by geometric predicates. The 
following predicate form for Simson’s Theorem is generated automatically 
from its constructive form mentioned above with MMP/Geometer. 

geom( [ [y5 ,x5 ,y4,x4,y3 ,x3 ,y2 ,x2 ,yl ,xl , v2 ,ul , vl] , [] , 

[A, [0,0] , B , [0,vl] ,C, [u2,v2] ,0, [xl,yl] ,D, [x2,y2] ,E, [x3,y3] , 

F, [x4,y4] ,G, [x5,y5]] , 

[ [cong , 0 , A , 0 , B] , [cong , 0 , A , 0 , C] , [cong , 0 , D , 0 , A] , 

[coll,E,A,B] , [perp,E,D,A,B] , [coll,F,B,C] , [perp,F,D,B,C] , 
[coll,G,A,C] , [perp,G,D,A,C]] , 

[ [sqdis , A,B] , [sqdis ,B , C] , [sqdis, A ,C] ] , [ [C0LL,E,F,G] ] ] ) ; 

Predicate [sqdis, A, B] means A ^ B with an algebraic representation as 
|OA| 2 = 0. The three predicates A / B,B yA C, A yA C are called the 
non-degenerated conditions (ndgs) for Simson’s theorem. 

A statement in predicate form can be represented by a 6-tuple: [mv,pv,pset, 
ps , ds,c\ where mv and pv are the main variables and parametric variables, 
pset is the set of points and their coordinates, ps is a set of predicates 
representing the hypotheses, ds is a set of predicates representing the non- 
degenerate conditions, c is the set of conclusions. 

Algebraic Form. In this form, coordinates are assigned to points in the state- 
ment and the hypotheses and conclusions are represented by algebraic equa- 
tions. It is straight forward to convert a statement in predicate form to 
algebraic form. The following is the algebraic form of Simson’s theorem in 
predicate form given above. 

geom( [ [y5 ,x5 ,y4,x4 ,y3 ,x3,y2 ,x2 ,yl ,xl , v2 ,ul , vl] , [] , 

[2*vl*yl-vl~2 , 2*v2*yl+2*u2*xl-v2~2-u2~2 , 
y2~2-2*yl*y2+x2~2-2*xl*x2,-vl*x3, 

-vl*y3+vl*y2 ,u2*y4-v2*x4+vl*x4-u2*vl , 
-v2*y4+vl*y4-u2*x4+v2*y2-vl*y2+u2*x2 ,u2*y5-v2*x5, 
-v2*y5-u2*x5+v2*y2+u2*x2] , 
[vl~2,v2~2-2*vl*v2+vl~2+u2'2,v2'2+u2'2] , 
[x4*y5-x3*y5-y4*x5+y3*x5+x3*y4-y3*x4] ] ) ; 
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Fig. 2. The relations between the four representation forms 



A statement in algebraic form can be represented by a 5-tuple: [ mv , pv, ps,ds, 
c] where mv and pv are the main variables and parametric variables, ps 
is a set of equations representing the hypotheses, ds is a set of equations 
representing the degenerate conditions, c is the set of conclusions. 

The relations between the four representation forms are illustrated in Fig. 2. 
An arrow from a form A to a form B means that MMP / Geometer can convert A 
to B. We use the method in [5] to convert a constructive form to a predicate 
form. The arrow marked with AGDG means that we will use AGDG methods 
reported in Section 4 to convert a predicate form to a constructive form. All 
other conversions are either obvious or can be found in [12]. 

A statement G = [ mv,pv,ps , ds, c] in algebraic form is said to be generically 
true if there exists a polynomial d in the variables pv such that 

Va; G mv, \/u G pv[(ps = 0Ads/0Ad^0)=tc = 0] 

If pv = 0, then G is universally true if 

Vc G mv[(ps = 0 A ds ^ 0) => c = 0] 

Since statements in all other forms can be converted into statements in al- 
gebraic form, we say a statement is generically true or universally true if its 
algebraic form is generically true or universally true. 

There are two ways to input a geometric statement graphically. 

Constructive. To input a statement constructively, we need to draw the dia- 
gram of the statement in a way similar to the ruler and compass construction. 
More precisely, we need to draw the objects in the statement sequentially 
with lines and circles. The drawing process can be naturally converted into 
a constructive description of the statement. For instance, to draw the dia- 
gram of Simson’s Theorem, we may just follow its constructive description 
given above to draw points A, B, C,0, D, E, F,G sequentially with functions 
provided by the software. 
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Declarative. In this form, we will first draw a sketch of the diagram and then 
add the geometric conditions. The software will adjust the diagram auto- 
matically such that these conditions will be satisfied. This process is called 
automated geometric diagram generation (AGDG). The declarative drawing 
process can be naturally turned into a predicate description of the statement. 
One way to draw a diagram in declarative form is to transform it into con- 
structive form automatically. Of course, not all statements in declarative 
form can be transformed into constructive form. For this kind of problems, 
we need to develop more complicated methods to draw their diagrams. This 
is the main task of AGDG, which will be introduced in Section 4. 

To draw the diagram of Simson’s Theorem in this way, we may first draw 
four points A, B , C , D on a circle and draw lines DE, DF, and DG such that 
E, F, G are on line AB , BC , AC respectively. Now we have a sketch of the 
diagram. Finally, we add the conditions DEEAB , DFEBC , DGEAC to 
the sketch and the software will adjust the sketch to satisfy these conditions. 

3 Automated Geometric Theorem Proving and 
Discovering 

3.1 Wu’s Method 

Wu’s method is a coordinate-based method for equational geometric statements, 
that is, the premise and the conclusion of the statement can be represented by 
algebraic or differential equations. The method first transfers geometric condi- 
tions into polynomial or differential equations, then deals with the polynomial 
equations with the characteristic set method [33, 36]. Variants of this method 
may be found in [2, 27, 30, 37]. 

Wu’s method is based on Wu-Ritt’s zero decomposition theorem for polyno- 
mial and differential polynomial equations [33, 36]. It may be used to represent 
the zero set of a polynomial equation system as the union of zero sets of equations 
in triangular form , that is, equation systems like 

fi(u,xx) = 0, f 2 (u,xi,x 2 ) = 0, . . . ,/ p (u, xi, ...,x p )= 0 

where the u could be considered as a set of parameters and the x are the variables 
to be determined. Let PS and DS be two sets of polynomials or differential 
polynomials, Zero (PS/DS) the set of solutions of PS = 0 over the field of 
complex numbers, which do not vanish any of the equations in DS. Wu-Ritt’s 
zero decomposition theorem [33, 30] may be stated as follows: 

Zevo(PS/DS) = U i Zero(SAT(A5 l )/D5) (1) 

where ASi are polynomial sets in triangular form and SAT(AS'i) are the satura- 
tion ideals of ASi. Variants of this algorithm may be found in [1, 2, 25, 26, 30, 37]. 
In MMP, the zero decomposition is implemented by Dingkang Wang [32]. 

Four versions of Wu’s method are implemented in MMP/Geometer. 
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WU-C. This method can be used to decide whether a geometric statement in 
constructive form is generically true , or it is true except some special cases called 
non-degenerate conditions (ndgs). The implementation is based on a variant of 
Wu’s method described in [5]. With this method, we may prove or dis-prove 
a geometric statement and give sufficient ndgs in geometric form. By sufficient 
ndgs, we mean 

(1) if the statement under consideration is valid in geometry textbooks (gener- 
ically true), then the statement is valid under these conditions; 

(2) if the statement is not valid under these conditions, then it will not become 
valid by adding more conditions unless the newly added conditions make the 
statement trivially valid or change the meaning of the statement. 

We use Simson’s Theorem as an example. 

wcproveC" Example Simson. Let D be a~point on the circumcircle 0 of 
triangle ABC. E is the foot from point D to line AB . F is the foot 
from point D to line BC. G is the foot from point D to line AC. 
Show that points E, F, and G are collinear . ") ; 

MMP/Geometer first converts the statement into constructive form , then proves 
the constructive statement and gives the following output. 

1. The statement is generically true, that is, it is true except some special cases. 

2. The sufficient ndgs are: — ■ [coll , A, B,C], A ^ B, B =£ C, A ^ C. 

WU-G. This method decides whether a geometric statement in algebraic form 
is valid [33, 36]. Let PS and DS be the equation part and inequation part 
corresponding to the ndgs of a statement. MMP/Geometer will first use Wu- 
Ritt’s zero decommission theorem to find triangular sets ASi as in (1). Let C = 0 
be the conclusion. If prem(C, AS/) = 0, then C = 0 is valid on Zero(SAT(ASi)). 
If ASi is irreducible, then this is also a necessary condition over the field of 
complex numbers. 

For Simson’s Theorem, we may prove that its following predicate form is 
valid. 

wprove ( [ [y5 ,x5 ,y4 ,x4,y3,x3 ,y2 ,x2,yl ,xl , v2 ,ul , vl] , [] , 

[A, [0,0] ,B, [0,vl] ,C, [u2,v2] ,0, [xl,yl] ,D, [x2,y2] ,E, [x3,y3] , 

F , [x4,y4] ,G, [x5,y5]] , 

[ [cong,0,A,0,B] , [cong,0,A,0,C] , [cong,0,D,0,A] , [coll,E,A,B] , 
[perp,E,D,A,B] , [coll,F,B,C] , [perp,F,D,B,C] , [coll,G,A,C] , 
[perp,G,D,A,C]] , [[sqdis,A,B] , [sqdis,B,C] , [sqdis,A,C]] , 

[ [COLL , E , F , G] ] ] ) ; 

Note that the result obtained here is stronger than that obtained with method 
WU-C: one ndg condition ->[coll, A, B 1 C] is removed from the description. 

WU-D. An advantage of Wu’s method is that it can be used to prove differential 
geometric theorems and mechanics [35]. The following is an example. 
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Example Kepler-Newton. Prove Newton’s gravitational laws using Kepler’s 
laws. Kepler’s first and second laws can be described as follows. 

Kl. Each planet describes an ellipse with the sun in one focus. 

K2. The radius vector drawn from the sun to a planet sweeps out equal areas 
in equal times. 

These laws can be expressed as the following differential equations. 

h x = r 2 - x 2 - y 2 = 0 

h 2 = a 2 - x" 2 - y " 2 = 0 

ki = r- p-ra = 0Ap' = 0Ae' = 0 

k 2 = x"y — y"x = 0. 



Newton’s law can be expressed as n\ = (ar 2 )' = 0. Then the problem is to prove 

(hi = 0 A h 2 = 0 A ki = 0 A k 2 = 0) => ni = 0. 

With MMP/Geometer, it is proved that the above statement is false. When add 
a ndg condition p ^ 0 (the ellipse does not becomes a straightline), we may use 
MMP/Geometer to prove the statement. 

depend ( [a , r , y , x] , [t]); 

wdprove ( [[a,r,y,x,p,e] ,[],[], [r~2-x~2-y~2 , a~2-x [2] ~2-y [2] "2, 
x*y [2] -x [2] *y , r-p-e*x] , [p] , [diff (a*r~2,t)]] ) ; 

Command depend([a,r,y, x\, [t]) defines a,r,y,x as functions in t. 

The following command proves that “A space curve C satisfies t = k! = 0 is 
a circle.” 

curve () ; wprove_curve ( [[],[],[], [t ,dif f (k,s)] , [] , 

[ [FIX_PLANE , C] , [FIX_SPHERE , C] ] ] ) ; 

The curve C , its arc length s, its torsion t and curvature k are defined in com- 
mand curveQ. diff(fc,s) is the differentiation of k with s. There are two con- 
clusions: [FIX_PLANE,C] meaning that C is on a plane and [FIX.SPHEREjC] 
meaning that C is on a sphere. 

WU-F. As pointed out by Wu [34], Wu-Ritt’s zero decomposition theorem can 
be used to discover geometric relations automatically. We use the following ex- 
ample to illustrate this. 

Example Heron-Qin Formula. Find the formula for the area of a tri- 
angle ABC in terms of its three sides. This problem may be solved with 
MMP/Geometer as follows 



wderive ( [[x,y,k] , [a,b,c] , [B, [0,0] , C,[a,0], A, [x,y]], 

[[dis, A, C, b] , [dis, A, B, c] , [area, k, A, B, C] ] , [] , [] ] ) ; 
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The above command will find the following relation between the last main vari- 
able k and the parameters: a, b , c automatically: 

16 * k 2 + a 4 — 2 * b 2 * a 2 — 2 * c 2 * a 2 + b 4 — 2 * c 2 * b 2 + c 4 = 0. 

In general, we may formulate the problem as follows. Let U = {u\, • • • , u m } 
be a set of parameters and X = {aq, • • • , x n } a set of dependent variables. The 
relation between the rq and the Xj is given by a set of algebraic or differential 
equations and inequations: 

P 1 (U,X) = 0,---,P s (U,X) =0 
Dr(U,X) ^0,---,D r (U,X) ^ 0. 



The problem is to find the relation between x\ and U, which can be done with 
Wu-Ritt zero decomposition under the following order u\ < U 2 ■ • • < u m < x±. 



3.2 The Area Method 



This method uses high-level geometric invariants such as the area and the 
Pythagorean difference as the basic tool of proving geometric theorems [8]. 
Instead of eliminating variables as in Wu’s method, the area method elimi- 
nates points from geometric invariants directly. The advantage is that short 
and human-readable proofs for geometric statements could be produced. This 
method works for constructive geometric statements. 

The signed area Sabc of triangle ABC is the usual area with a sign depending 
on the order of the three vertices of the triangle. The following properties of the 
signed area is clearly true. 



A1 Let M be the intersection of two non-parallel lines AB and PQ and Q ^ M. 
Then 2M = |£^.. 

QM OQAB 

A2 If Aj^B, PQ || AB iff Spab = Sqab- 



Let us consider the following example. 



aprove (" Example Parallelogram. Let ABCD be a parallelogram. 0 is 
the intersection of diagonals AC and BD . Show that 0 is the 
midpoint of AC."); 

MMP/Geometer produces the following proof for the parallelogram theorem, 
which is with clear geometric meaning. The comments on the right hand side is 
added by the authors. 



AO Sabd 

OC Sbcd 

S a no -i 

Sabc 



By Al. 

By A2, Sabd = Sabc , Sbcd = Sabc- 



To deal with perpendicularity, we need to define the Pythagorean difference. 
For points A, P, and C, the Pythagorean difference, Pabc , is defined to be 
Pabc = AB + CB — AC . The following Pythagorean theorem is taken as 
a basic (unproven) property of the Pythagorean difference. 
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Fig. 3. Parallelogram Theorem and Orthocenter Theorem 



PI AB±PQ iff P PAB = Pq AB - 

Let us consider the following example. 

aproveO Example Orthocenter. ABC is a~triangle. F is the foot from 
A~to BC. E is the foot from B to AC. H is the intersection of 
lines AF and BE. Show that CH is perpendicular to AB."); 

By PI, we need only to show Pach = Pbch- MMP/Geometer produces the 
following proof for the Orthocenter theorem. 

Pa oh Pacb Paoh __ ^ 

Pbch Pbca Pacb 

By PI, we can eliminate point H as follows: Pach — Pacb ; Pbch = Pbca- 



The following is the machine proof for Simson’s Theorem. Point G\ = EF n 
AC is added by MMP/Geometer automatically. We need only to prove — = 

C Ct 

==■ Details on this proof may be found in [8] 

CGi 



The machine proof 



(4£)/(AG^) 
v BG 7 ' v BGi 7 



^ SBEF . AG 
S AEF BG 

EL PbAD'SbEF 
Saef-(-Pabd) 



E —Pbad-Pacd-SabePaca 
(-Pcad-Sace)-Pabd-Paca 

simplify P BAn .P Ann -S ARK 
P CAD S ACE- p ABD 

E Pbad-Pacd-Pcbd-Sabc-Pbcb 
Pcad(-PbcdSabc)Pabd-Pbcb 



simplify P RAn .P Ann .P nFin 
-Pcad-Pbcd-Pabd 



The eliminants 

AG-j 'WL S A ef — £ Praf> 

BGi S BEF 5 BG ~ P ABD 



S AEF = 



F -Pcad-Sace 



P AC A 



S BEF = 



F Pacd -Sabe 



P AC A 



S ACE 



E -PrOD-SaBC 
P BCB ’ 



S ABE 



E P C BD S ABC 
P BCB 



P ABD= - 2 (bd-ab ■COS (AD)) 
Pbcd= — 2(CDBC-cos(BD)), 

p CAD=^)aD-AC •COS (CD)) 

p C bd=^{bdbc •COS (CD)) 
P ACD= - 2 {CD AC •COS (AD)) 
p bad=2(adab ■COS (BD)) 



co—cir (2 AD-AB-cos(BD))(-2CD-AC-cos(AD))-(2BDBC-cos(CD)) 
-(2ADAC-cos(CD))-(-2CDBC-cos(BD))-(-2BDABcos(AD)) 



simplify 



1 
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A E B 

Fig. 4. Orthocenter Theorem 



3.3 Deductive Database Method 

For a given geometric diagram, MMP/Geometer can generate a geometric de- 
ductive database (GDD) which contains all the properties of this diagram that 
can be deduced from a fixed set of geometric axioms, and for each geometric 
property in the database, MMP/Geometer can generate a proof in traditional 
style [9]. 

Let Do be the premise of a geometric statement and R the set of geometric 
axioms or rules. We may use the breadth-first forward chaining method to find 
new properties of the corresponding diagram. Basically speaking, the method 
works as follows 




R 

C 




R 

... c 



Dk 



(Fixpoint) 



where |_i is the union of Di and the set of new properties obtained by applying 
rules in R to properties in Di. If at certain step Dk = Dk+ 1 , i.e. , 



R{D k ) = D k , 

then we say that a fixpoint (of reasoning) for Do and R is reached. 

The naive form of breadth-first forward chaining is notorious for its ineffi- 
ciency. But, in the case of geometric reasoning, by introducing new data structure 
and search techniques, we manage to build a very effective prover based on this 
idea [9]. 

The GDD method may take a statement in predicate, constructive, or natural 
language form as input. Let us consider the following form of the Orthocenter 
Theorem (Fig. 4). 

prove_gdd( " Example Orthocenter . ABC is a~triangle. D is the foot 
from B to AC. E is the foot from C to AB. F is the intersection of 
lines BD and CE. G is the intersection of lines BC and AF."); 

The initial database is the hypotheses: D, A, C ; E, A, B\ F, B, D; F, C , E\ G, 
B,C ; G,A,F are collinear sets, BD _L AC,CE _L AB. The fixpoint contains 
151 geometric properties: 

6 collinear point sets, 
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3 perpendicular pairs, 

6 co-cyclic points sets, 

24 equal angle pairs, 

1 similar triangles sets, 

105 and equal ratio pairs. 

The forward chaining is a natural way of discovering properties for a given 
geometric configuration. Any thing obtained in the forward chaining may be 
looked as a “new” result. Take the simple configuration (Fig. 4) related to the 
orthocenter theorem as an example. MMP/Geometer has discovered the most 
often mentioned properties about this configuration: AG _L BC (the orthocenter 
theorem) and LEG A = LAGD. The fixpoint also contains seven sets of similar 
triangles 



A DBA - 


- A DCF r- 


- A EBF - 


- A EC A] 


A DCB r 


- A DFA r 


- A GFB - 


- A GCA- 


AEFA - 


- A EBC r. 


- A GBA - 


* A GFC-, 


A FBC ~ 


- A FED - 


- A GBD r 


« AGEC; 


AACB - 


- A AED - 


- A CCD r. 


- A CEB- 


A CED r 


- A CAF - 


- A GAD r. 


- A GEF; 


A FBA - 


- A EBD r. 


- A FGD - 


- A EGA. 



For each geometric property in the database, MMP/Geometer can produce a 
proof of traditional style. The following is the proof for the Orthocenter Theorem 
(. AG T BC) by MMP/Geometer. Notice that the proof is in an “analysis style”, 

i.e. , it starts from the conclusion and goes all the way to the hypotheses of the 
statement. 

1. AG _L BC, 

because AC T BD (hypothesis), (2) L [AC, BD] = L[BC,AF\. 

2. L[AC,BD\ = L[BC,AF ], 

because (3 )L[AC, BC] = L[BD,AF ].( This is rule in MMP/Geometer ). 

3. L[CA,CB } = L[BD,AF), 

because (4 )L[CA,CB] = L[DE,AB], (5 )L[BD,AF] = L[DE,AB}. 

4. L[CA,CB } = L[DE, AB], 

because (6) co-cyclic [.B, D, C, E], 

5. L[BD, AF] = L[DE, AB\, 

because (7) co-cyclic [A, D, E, F], 

6. co-cyclic[B,D,C,E], 

because DC T DB (hypothesis), EC _L EB (hypothesis). 

7. co-cyclic[A, D, E, F], 

because DF T DA(hypothesis), EF T EA(hypothesis). 
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The first step of the proof can be understood as follows. AG _L BC is true 
because AC _L BD which is a hypothesis and L [AC, BD] = /.[BC, AF] which 
will be proved in the second step. The other steps can be understood similarly. 

One of the largest databases obtained with MMP/Geometer is for the nine- 
point circle theorem (Fig. 5(a)), which contains 6019 geometric relations. In 
predicate form, the database contains 6646428 predicates. 

prove_gdd(" Example Nine-point-circle. ABC is a triangle. K is the 
midpoint of BC. L is the midpoint of CA. M is the midpoint of 
AB . D is the foot from B to AC. E is the foot from C to AB. F is 
the intersection of lines BD and CE. G is the intersection of 
lines BC and AF. H is the midpoint of AF. I is the midpoint of 
BF. J is the midpoint of CF.Show that K, G, J, D are 
cocircle . ") ; 

Constructing new points or lines is one of the basic methods for solving geo- 
metric problems. One advantage of MMP/Geometer is that it can automatically 
add auxiliary points to prove a geometric statement if needed [9] . The following 
example is shown in Fig. 5(b). 

prove_gdd( "Example Trapezoid. Let ABCD be a trapezoid. M is the 
midpoint of AC. N is the midpoint of BD. E is the intersection of 
lines MN and BC. Show that E is the midpoint of BC."); 

The conclusion is not in the first fixpoint. The program then automatically 
adds an auxiliary point A$ which is the midpoint of AD. With this auxiliary 
point, MMP/Geometer generates a fixpoint containing the conclusion. 

The Butterfly theorem in Fig. 5(c) also needs auxiliary points. 

prove_gdd( " Example Butterfly-Theorem. P, Q, R, and S are four 
points on a circle 0. A is the intersection of lines PQ and SR. N 
is the intersection of line PR and the line passing through A and 
perpendicular to 0A. M is the intersection of line QS and the 
line passing through A and perpendicular to 0A. Show that A is 
the midpoint of NM."); 
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Wu’s Method 


Area Method 


GDD Method 


Example 


Time 


Example 


Time 


Example 


Time 


Simson (WU-C) 


20 


Parallelogram 


10 


Orthocenter 


110 


Simson (WU-G) 


80 


Orthocenter 


10 


Trapezoid 


50 


Kepler-N ewton 


3200 


Simson 


30 


Butterfly 


110 


Curve 


2100 






Nine-point 


560(s) 


Herron-Qin 


10 











Table 1. Statistics for the examples 



For this problem, MMP/Geometer automatically adds an auxiliary point Ao 
which is the intersection of the line passing through S and parallel to AN and the 
circle O. With this point, MMP/Geometer generates a fixpoint which contains 
the conclusion. 

3.4 Implementation and Experimental Results 

MMP is developed with VC in a Windows environment. Current version may be 
found in the webpage listed in [12]. Many techniques are introduced to enhance 
the efficiency of the proving algorithms. One of the main operation in Wu’s 
method is to test whether the successive pseudo-remainder of a polynomial P 
with respect to a triangular set TS is zero. We implements three techniques to 
enhance the efficiency. 

Simplifying. Let TS be the triangular set and LT S the set of the polynomials 
in TS which contain less than three terms and are linear in their leading 
variables. Then the pseudo remainder of a polynomial Q with LT S contains 
equal or less terms than Q. Therefore, when we take the pseudo remainders 
of P and the polynomials in T S with respect to LT S to obtain a polynomial 
P' and a new triangular set TS and then compute the pseudo-remainder 
of P' with respect to TS'. 

Branching. Suppose that we need to compute the pseudo-remainder prem(c, d) . 
If c = c s x s k + .. + Co, d = dtx * + .. + do and j < k , then we can check whether 
prem(ci, d ) = 0,i = 0,...,s instead of whether prem(c,d) = 0. Generally 
speaking, each prem(ci,d) has fewer terms than prem(c, d) and is easy to 
compute. 

Dividing Extraneous Factors. It is known that extraneous factors appear in 
the process of doing pseudo remainders. These factors are removed automat- 
ically to reduce the size of the polynomials in the proof. 

Table 1 contains the running times for the examples in this Section 3. The 
data is collected on a PC compatible with 1.5G CPU. The time is in microsec- 
onds. 

Why do we implement more than one method in MMP/Geometer? First, each 
method has its advantages and shortcomings. Generally speaking, the proving 
power of the methods are as follows 
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Methods 


No. of Theorems 


Source 


WU-C 


512 


[2] 


WU-G 


450 


[3] 


WU-D 


100 


[6] 


WU-F 


120 


[4] 


AREA 


400 


[8] 


DBASE 


170 


[9] 



Table 2. Theorems proved using Wu’s method, area method and GDD method 



WU-C > WU-G > AREA > DBASE. 

Table 2 gives the number of theorems proved with these methods. When consid- 
ering to produce elegant and human-understandable proofs, the order is reversed. 

Second, with these methods, for the same theorem, the prover can pro- 
duce a variety of proofs with different styles. This might be important in using 
MMP/Geometer to geometry education, since different methods allow students 
to explore different and better proofs. 

Only a selective set of examples shown in Table 2 was tested in 
MMP/Geometer. We expect that almost all of the theorems mentioned in the 
above table can be proved with MMP/Geometer since the methods used in 
MMP/Geometer are improved version for the original methods used to prove 
them. 



4 Automated Geometric Diagram Generation 

4.1 Dynamic Geometry 

By dynamic geometry , we mean computer generated geometric figures which 
could be changed dynamically. Generally, we may perform the following opera- 
tions on the figures: dynamic transformation, dynamic measurement, free drag- 
ging, and animation. By doing dynamic transformations and free dragging, we 
can obtain various forms of diagrams easily and see the changing process vividly. 
Through animation, the user may observe the generation process for figures of 
functions. 

Most dynamic geometry software systems [23, 22, 13] use construction se- 
quences of lines and circles to generate diagrams. Since such sequences are easy 
to compute, dynamic geometry software systems are usually very fast. 

As an example, let H be the orthocenter of triangle ABC. We fix points A 
and B and let point C move on a circle c. We want to know the shape of the locus 
of point H. Let (xo,yo) and r be the center and radius of circle c, A = (0,0), 
and B = (d, 0). With Wu’s method of mechanical formula derivation, we may 
obtain the equation of this locus 



((x - x 0 ) 2 + Vq- r 2 )y 2 + 2y 0 (x - d)xy + (x - d) 2 x 2 = 0 
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with MMP /Geometer as follows. 

wderiveC [[y3,x3,y2,x2,yl,xl,y] , [x,d,x0,y0,r] , 

[A, [0,0] ,B, [d,0] ,0, [x0,y0] ,C, [xl,yl] ,H, [x,y]] , 

[ [dis , C , 0 , r] , [perp , A , H , B , C] , [perp , B , H , A , C] ],[],[]]) ; 

But from this equation, we still do not know the shape of the curve. With 
MMP/Geometer, we can draw the diagram of this curve as the locus of point H. 
By continuously changing the radius of circle c, we may observe the shape 
changes of the curve (Fig. 6). 

The diagram in Fig. 7 is to generate the locus of the moon when the moon 
rotates around the earth on a circle and the earth rotates around the sun on an 
ellipse. For its description, please consult [12]. 

In dynamic geometry, the drawing process is based on the constructive de- 
scription for a diagram, that is, each geometric object is constructed with ruler 
and compass. In a construction sequence, free and semi-free points could be 
dragged. The computation can be carried out in real-time because the equations 
raised from a construction sequence are almost in triangular form: variables are 
introduced at most two by two. Further, only equations with degree less than 
or equal to two are involved. MMP/Geometer has all the above functions of 
dynamic geometry. 




Fig. 6. Continuous change for the locus of the orthocenter 




Fig. 7. Locus of the moon 
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A diagram with 21 given distances 



4.2 Intelligent Dynamic Geometry 

Most geometric diagrams in geometry textbooks are described declaratively, 
and the task of converting such a description to constructive form is usually 
done by human. For some diagrams, it is difficult to find a constructive solu- 
tion, and for more diagrams there exist no ruler and compass constructions. In 
MMP/Geometer, by combining the idea of dynamic geometry and AGDG we 
implement an intelligent dynamic geometry software system, which can be used 
to input and manipulate diagrams conveniently. 

There are two major steps to draw a declaratively given diagrams. First, 
we try to find a ruler and compass construction by mechanizing some of the 
techniques of ruler and compass construction developed since the time of ancient 
Greek. If we fail to do so, general AGDG methods are used to draw the diagram. 

To find a ruler and compass construction for a geometric diagram in predicate 
form, we first transform the geometric relations into a graph and then solve it 
in two steps. 

1. We repeatedly remove those geometric objects that can be constructed ex- 
plicitly until nothing can be done. This simple algorithm is linear [11] and 
solves about eighty percent of the 512 problems reported in [2]. 

2. If the above step fails, we use Owen and Hoffmann’s triangle decomposition 
algorithms [20, 29] to reduce the problem into the solving of triangles. 

3. If the above step fails, we use certain geometric transformations to solve 
the problem[Ll]. This is a quadratic algorithm and is complete for drawing 
problems of simple polygons. 

To see the power of the methods, let us look at the example in Fig. 8, where 
each line segment represents a given distance between its two end points. We 
have twenty one distances. In order to solve the problems algebraically, we need 
to solve an equation system consisting of twenty one quadratic equations. But 
with the triangle decomposition [20, 29], the problem can be reduced to solving 
of triangles, as shown in Figures 8(b) and (c). As a consequence, the problem is 
ruler and compass constructible. 

As another example, let us draw the quadrilateral in Fig. 9, which cannot be 
solved with the triangle decomposition method. This diagram cannot be drawn 



62 



Xiao-Shan Gao and Qiang Lin 



^L-_ 



p 




P2^ 




Pj ^ 6 ^Ps 



Fig. 9. Lengths of four edges and angle (L 2 ,L4) are given 



with ruler and compass explicitly. To solve it, MMP/Geometer adds a parallelo- 
gram P1P2P3P2 [ 11 ]. Since ZP4P1P2, l-Pi-P^li | P'i P'i are known, it is easy to find 
a construction sequence for quadrilateral P3P4P1P2. Point P 2 can be constructed 
easily. 

Using methods of AGDG allows us to have more power to manipulation 
the diagram. If a construction sequence for a diagram has been given, we may 
only drag the free and semi-free points in the diagram. This drawback may be 
overcomed as follows. Suppose that we want to drag a point, we may re-generate 
a new construction sequence in which the point is a free point. This kind of 
dragging is called intelligent dragging. 

Let us consider the diagram in Fig. 4 . A construction order for the points in 
this diagram is as follows 



where A, B , C are free points. F is the intersection of lines BD and CE , hence 
a fixed point. In dynamic geometry software, we cannot drag this point. But 
in MMP/Geometer, when a user wants to drag this point, a new construction 
sequence 



is automatically generated, in which F is a free point and can be dragged. 

If we cannot find a ruler and compass construction sequence for a problem, 
MMP/Geometer will try to solve the problem with the following algorithms. 

1 . Graph analysis methods are used to decompose the problem into c-trees [ 14 ]: 

2 . The c-tree is reduced to a general construction sequence 



where C/ are sets of geometric objects such that 
— Ci can be constructed from UjjlyCfc. 

— and Ci is the smallest set satisfying the above condition. 

3 . Compute the position of Ci from Ci, ■ ■ • , Cj_i. MMP/Geometer uses two 
methods to do this. 

Optimization Method The problem is converted to solving a set of alge- 
braic equations: 






/i(X) = 0, • • • , / m (X) = 0 
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Fig. 10. Packing circles into a triangle 



where X is a set of variables. Let 

m 

<r(x) = £/?. 

i — 1 

We use the BFGS method to find a position Xo such that er(Xo) is 
a minimal value [17]. If cr(Xo) = 0 then Xo is a set of solutions for 
original equation system. 

Locus Intersection Method (LIMd) The above method based on opti- 
mization can find one solution only. If we want to find all the solutions, 
we may use the LIMd method [15]. The LIMd method is a hybrid method 
to find all solutions to geometric problems. 

Let us consider a problem of packing circles into a fixed triangle. The prob- 
lem is to pack n(n + l)/2 circles (n rows of circles) tangent to adjacent circles 
and the adjacent neighboring sides of a given triangle. Fig. 10 is the cases for 
n = 1, 2, 8. The most difficult case is n = 8, in which we need to solve 24 quar- 
tic equations and 84 quadratic equations simultaneously. This equation system 
cannot be simplified essentially. Table 3 shows the running times for different n 
of this problem on a PC with CPU 1.5G. It seems that problems of large size 
can be solved in real time. 



# circles (# rows) 


# quartic equations 


# quadratic equations 


time 


3 (2) 


6 


3 


0.18 


6 (3) 


9 


9 


0.50 


10 (4) 


12 


18 


1.72 


15 (5) 


15 


30 


1.93 


21 (6) 


18 


45 


4.15 


28 (7) 


21 


63 


10.92 


36 (8) 


24 


84 


15.45 



Table 3. Running statistics for the tangent packing problems 



64 



Xiao-Shan Gao and Qiang Lin 



4.3 Implementation Issues 

After an algorithm is selected, we still need to add details in the implementation 
to enhance the performance of the algorithm. One problem is the ambiguities 
rising from operations like “intersection of two circles” or “intersection of a circle 
and a line”. These constructions have several solutions. When the user drags 
a point randomly, the program will compute the position of all points in the 
figure automatically. Our goal is to keep a continuous movement of dependent 
elements. In other words, we try to avoid “jumps”. The problems is solved by 
comparing the two solutions with the initial positions and then the software will 
remember the relative position of the relevant elements. 

During the generation of construction sequences or generalized construction 
sequences, we need to select a starting elements in order to obtain better con- 
struction sequences. We adopt the following strategies: (1) Use the input order 
by the user as the initial order; (2) always try to construction lines before points; 
(3) try to use the distance constraints before other constraints. 

Another major concern is the speed. While we dragging the diagram ran- 
domly, we hope to see smooth movements. We adopt the following strategies to 
keep the software fast: (1) if the dragging element o is not the first element in 
a construction sequence, then we need only to compute the positions for those 
elements after o in the construction sequence; (2) when searching construction 
sequences, try to use the constructions which are easier to compute; (3) opti- 
mize the code by using explicit formulas to solve linear and quadratic equations. 
The result is encouraging: the software may give smooth movements for most 
problems containing less than one hundred geometric elements. 

5 Conclusions 

MMP/Geometer is a geometric software. It can be used to prove and discover 
theorems in Euclidean and differential geometries. It can also be used to pro- 
duce proofs with geometric meanings and proofs in traditional style. We are 
currently building a webpage which will include MMP/Geometer and a large 
set of geometric theorems including those mentioned in Table 2. 

MMP/Geometer can also be used as a geometric diagram editor. For a given 
input, MMP/Geometer will draw the diagram first using ruler and compass con- 
struction and then using general AGDG methods. It seems that large diagrams 
such as those in Fig. 10 can be handled satisfactorily. 

Altogether, we hope that MMP/Geometer may provide a useful tool for 
people to study, to learn and to use geometry. 
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Abstract. Formalized proof schemes are the starting point for testing, 
comparing, and benchmarking of different geometry theorem proving ap- 
proaches and provers. To automatize such tests it is desirable to collect 
a common data base of proof schemes, and to develop tools to extract 
examples, prepare them for input to different provers, and run them 
“in bulk”. The main drawback so far of special collections, e.g., Chou’s 
collection [2] with more than 500 examples of proof schemes, was their 
restricted availability and interoperability. We report about first experi- 
ence with a generic proof schemes language, the GeoCode language, that 
was invented to store more than 300 proof schemes in a publicly available 
repository, and tools to prepare these generic proof schemes for input to 
different target provers. 

The work is part of the SymbolicData project [15]. 



1 Introduction 

Proofs in classical textbooks on geometry are usually tricky, involve a big portion 
of heuristics and are hard to formalize and even to systematize. Nevertheless 
just these circumstances - the small set of prerequisites necessary to formulate 
the problems and the erudition and non formal approaches required to solve 
them - generate the great interest even of leisure mathematicians in geometry 
problem solving. With this experience in mind it seems to be impossible to find 
a systematic way to prove geometry theorems. 

But such a unified approach provides, e.g., the coordinate method: Translate 
geometric configurations into algebraic relations between coordinates and try to 
solve the algebraic counterpart of the geometric problem by algebraic methods. 
It was this framework that inspired the young Gauss for his famous solution to 
construct a regular 17-gon by ruler and compass. 

With increasing capabilities of modern computer equipment to perform sym- 
bolic algebraic manipulations the attempts to algorithmize and mechanize this 
part of mathematics made great progress. Many deep results on mechanized ge- 
ometry theorem proving were discovered by the “Chinese provers” in the school 
of W.-T. Wu at MMRC and their collaborators all over the world. For a survey 
on the current state of the art see, e.g., the monographs [2, 3, 4, 21, 22] or the 
conference proceedings on “Automated Deduction in Geometry” [ .7, 5, 14] and 
this volume. 



F. Winkler (Ed.): ADG 2002, LNAI 2930, pp. 67-86, 2004. 
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The proofs are not “automated” but “mechanized” in the following sense: 
A partly informal human readable geometric statement requires a translation 
into proof schemes or coordinatizations written in a strong computer readable 
syntax. These (yet purely geometric) proof schemes are the starting point for 
automated application of different proving approaches. Although straightforward 
in most cases, some theorems require a quite tricky coordinatization to get an 
algebraic translation that can be handled by current provers and computers. 

To test, compare and benchmark different provers and software it is desirable 
to have such proof schemes collected in a common publicly available repository 
in digital form and a common format that easily translates to the special input 
formats of different provers. This requires to agree upon a generic proof scheme 
language standard that can be managed with appropriate tools by all interested 
parties. 

Below we describe a first approach to such a standard, the GeoCode. It is 
used to store proof schemes as GEO records in a repository that is publicly 
available as part of the SymbolicData project [15]. At the moment the Symbolic- 
Data GEO collection contains more than 250 proof schemes, mainly from [2]. 
Special Perl based tools support the syntactical translation of GeoCode into 
the special input format of provers that implemented this common interface. At 
the moment - as a first reference application - this interface is implemented in 
the different versions of the author’s GeoProver packages [8], that provide tools 
to run proof scheme translations based on the coordinate method on one of the 
major CAS (Maple, Mathematica, MuPAD, Reduce). Figure 1 shows the relation 
between the SymbolicData PROBLEMS and GEO records, the GeoCode syntax 
definition, the GeoProver (lower row), and different stages of the process of 
mechanized theorem proving (upper row). 

This paper starts with some background on geometry theorem proving (sec- 
tion 2). Then we describe the design of the GEO records, the syntax of the 
GeoCode standard and the GeoProver packages as an implementation of that 
standard (section 3). In section 4 we discuss by means of examples how to com- 
pile new (generic) proof schemes, to translate them into GeoProver notion, to 
run this code on different CAS and to experiment with the resulting algebraic 
problems. 

For real usability of the GeoCode concept one has to estimate the efforts 
required to implement this standard for other provers. There is not yet practi- 
cal experience but the semantic similarity of “foreign” special geometric proof 
schemes suggests that such an interface should easily be implemented. The prob- 
lem is discussed to some extend in section 5. 

2 The Coordinate Method 

The main approach to mechanized geometry theorem proving considered so far 
uses the coordinate method: It translates geometric statements into their al- 
gebraic counterparts, i.e. , statements about systems of polynomial or rational 
functions, and tries to solve these algebraic problems by algebraic methods. 
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Fig. 1. Formal and informal elements in mechanized geometry theorem proving 
and their counterparts in the SymbolicData concept 

A particularly simple proof technique applies to geometric configurations of 
constructive type. To prove, e.g., the median intersection theorem , take three 
generic points A, B , C, construct the midpoints Ai,Bi,Ci of the segments BC, 
AC and AB and check for the lines AA \ , BB\ and CC\ to satisfy the concurrency 
condition. Here is this proof scheme in GeoCode syntax: 

<Points> 

$A : =Point [al , a2] ; $B : =Point [bl ,b2] ; $C : =Point [cl , c2] ; 

</Points> 

<coordinates> 

$A1 : =midpoint [$B, $C] ; $B1 : =midpoint [$A, $C] ; $C1 : =midpoint [$A, $B] ; 
</coordinates> 

<conclusion> 

is_concurrent [pp_line [$A, $A1] , pp_line [$B , $B1] , pp_line [$C , $C1] ] ; 
</conclusion> 

or even shorter 

<Points> 

$A : =Point [al , a2] ; $B : =Point [bl ,b2] ; $C : =Point [cl , c2] ; 

</Points> 

<conclusion> 

is_concurrent [median [$A,$B,$C] , 

median [$B , $C , $A] , median [$C , $A , $B] ] ; 



</conclusion> 
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This proof scheme has pure geometric nature and can serve also for input to 
provers not based on the coordinate method or, to some extend, even for diagram 
drawing tools. With SymbolicData Perl tools the proof scheme translates, e.g., 
to the GeoProver MuPAD syntax 

A : =Point (al , a2) ; B : =Point (bl ,b2) ; C : =Point (cl , c2) ; 

A1 : =midpoint (B ,C) ; Bl : =midpoint (A,C) ; Cl : =midpoint (A,B) ; 
is_concurrent(pp_line(A,Al) , pp_line(B,Bl) , pp_line(C,Cl)) ; 

where now (the GeoPrcwer-implementation of) midpoint(X, Y) computes the 
coordinates of the midpoint of the segment XY and pp_line(X, Y) those of 
the line through X and Y. is_concurrent(a, 6, c) translates to a polynomial 
condition in the coordinates of a,b,c (in fact, a determinantal expression) that 
vanishes iff these lines pass through a common point. The return values of all 
these functions are (tuples of) rational expressions in the coordinates of the 
formal input parameters. To prove this theorem (and other theorems of this 
kind) means to compose a nested rational expression and to check if it simplifies 
to zero. 

In general, we say that a geometric configuration is of constructive type , if 
its generic configuration can be constructed step by step in such a way, that 
the coordinates of each successive geometric object can be expressed as ratio- 
nal functions in the coordinates of already constructed objects and algebraically 
independent variables, and the conclusion can be expressed as vanishing of a ra- 
tional function in these coordinates. Note that due to Euclidean symmetry even 
for generic configurations some of the coordinates may be chosen in a special 
way. 

Geometry theorems with non-linear geometric conditions that cannot be 
solved rationally in the given indeterminates require other proof techniques. 
Consider, e.g., the angular bisector intersection theorem: 

The condition on_bisector(P,A,B,C) for a point P = Point(xi, X 2 ) to be 
on either the inner or the outer bisector 1 of L ABC translates into a polynomial 
of (total) degree 4 in the generic coordinates of A, B , C and quadratic in the 
“dependent” coordinates of P. 

To prove the angular bisector intersection theorem we collect the conditions 
on P to be on each of the bisector pairs through A and B, solve the corresponding 
polynomial system for the coordinates of the intersection point (s) P and show 
that each solution belongs also to (one of) the bisectors through C. 

In more detail: Due to Euclidean symmetry we can choose special coordinates 
for A and B to simplify calculations. 

A : =Point (0 , 0) ; B : =Point (1 , 0) ; C : =Point (ul ,u2) ; P:=Point(xl,x2) ; 
polys:={on_bisector(P,A,B,C) , on_bisector(P,C,A,B)}; 
con:=on_bisector(P,B,C,A) ; 

1 There is no way to distinguish between the inner and outer bisectors in unordered 
or dynamic geometry. 
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polys is a system of two polynomial equations of degree 2 in (aq,^) with 
coefficients in Q(mi,u 2 ). It has 4 (generic) solutions that correspond to the 4 
(generic) intersection points of the bisector pairs through A and B. They can be 
computed, e.g., with Maple: 

sol : =solve (polys , {xl ,x2}) ; 



The solution involves algebraic RootOf -expressions that require a powerful 
algebraic engine to cope with. For Maple, simplify (subs (sol, con)) returns 0 
and thus proves the theorem generically. 

We can also reformulate the geometry theorem as a vanishing problem of the 
polynomial conclusion on the zero set of the system of polynomials that describe 
the given geometric configuration and solve it with a Grobner basis and normal 
form calculation (Maple) 

with(Groebner) : 

TO : =plex (xl , x2) : gb : =gbasis (polys , TO) : 
normalf (con,gb,T0) ; 

In general, this kind of algebraization of geometry theorems yields a polyno- 
mial ring S = fc[v] with variables v = (iq, . . . , v n ), a polynomial system F C S 
that describes algebraic dependency relations in the given geometric configura- 
tion, a subdivision v = x U u of the variables into dependent and independent 
ones, and the conclusion polynomial <?(x, u) £ S. To prove the corresponding ge- 
ometry theorem t/(x,u) must vanish on all “geometrically relevant” components 
of the set of common zeroes of F. We get proof schemes of equational type. 

If we compute in the ring So = fc(u)[x] as we did in the above example, 
exactly those components of the zero set of F are visible where the variable set 
u remains independent. Hence if the normal form of g wrt. a Grobner basis G 
of F computed in So vanishes, as in the example above, the geometry theorem is 
generically true , i.e., holds if some “special” configurations for the independent 
parameters are avoided that correspond to degenerate geometric situations. More 
subtle examples can be analyzed with the Grobner factorizer or more advanced 
techniques. 

There are also other algebraic techniques to analyze such polynomial systems, 
e.g., based on pseudo division and triangular sets. See [12] or the monographs [4, 
20] for a survey. 




%1 = RootOf (4 U 2 -if 4 + (—8 Mi 2 — 8 U 2 2 + 8 ui) -Z 3 

+ (—4 mi u 2 + 4 u\ 2 U 2 — 4 u 2 + 4 u 2 3 ) -Z 2 + 4 u 2 2 -Z — u 2 3 ) ■ 
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Formalized proof schemes are the starting point for testing, comparing, and 
benchmarking different geometry theorem proving approaches and provers. To 
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automatize such tests it requires a data base of proof schemes, and tools to 
extract examples, prepare them for input to provers under consideration, and 
run them “in bulk” . 

To compile such a data base a “proof writer” has to fix (realistic) proof 
schemes for given informal statements of geometry theorems. Several people 
collected such data, e.g., S.-C. Chou who discussed in [2] more than 500 examples 
of geometric statements and appropriate algebraic translations. 

It is desirable not do that work twice but compile proof schemes in such 
a way that they can be easily prepared for input to different provers. The main 
drawback so far of special collections is the restricted interoperability of proof 
schemes. To fix a proof scheme for automated processing by different provers 
requires a generic language that can be mapped to all target systems. Below 
we report on our experience with the GeoCode language that was developed to 
store generic proof schemes in the SymbolicData GEO record collection. 

During our work on that collection we stored (and partly modified and 
adapted) about 200 of them. We collected also solutions of geometry problems 
from other sources, e.g., the IMO contests, see [11]. Much of this work was done 
by my “proof writers” , the students Malte Witte and Ben Friedrich, who com- 
piled first electronic versions for many of these examples. 

3.1 The SymbolicData Project 

The SymbolicData project was set up to create and manage a publicly available 
repository of digital test and benchmark data from different areas of symbolic 
computation and to develop tools and concepts to manage such data both in the 
repository and at a local site. In a first stage we concentrated on the development 
of practical concepts for a convenient data exchange format, the collection of 
existing benchmark data from two main areas, polynomial system solving and 
geometry theorem proving, and the development of appropriate tools to process 
this data. A tight interplay between conceptual work, data collection, and tools 
(re)engineering allowed continuously to evaluate the usefulness of each of the 
components. 

For easy reuse we concentrated on free software tools and concepts. The data 
is stored - one record per file - as tag/value pairs in a XML like ASCII format 
that can be edited with your favorite text editor. The tools are completely written 
in Perl using modular technology. The data format definitions are part of the 
data base itself and stored in a similar way. It can be specified and extended by 
the user in an easy manner and very broad range to add new material and/or 
to modify existing one. I refer to [1, 9, 10] and the SymbolicData documentation 
for more details. 

The project is organized as a free software project with a CVS repository 
equally open to people joining the SymbolicData project group. The Symbolic- 
Data project is part of the benchmark activities of the German “Fachgruppe 
Computeralgebra” who also sponsored the web site [15] as a host for presentation 
and download of the tools and data developed and collected so far. We kindly 
acknowledge support also from UMS MEDICIS of CNR/Ecole Polytechnique 
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(France) who provides us with the needed hard- and software to run this web 
site. 

Tools and data are available for local installation (without any registration) 
under the terms of the GNU Public License as tar-files from our Web site. 

3.2 The Structure of the SymbolicData Proof Scheme Records 

Each SymbolicData record has some common attributes (Type. Key, CRef , . . . ) 
that serve for identification or store relational information. The other attributes 
store different parts of the proof scheme in GeoCode syntax and will be described 
now. See [10] for a detailed description of the overall GEO record structure. 

The SymbolicData GEO proof scheme records are divided (roughly) into two 
types according to their prooftype attribute: constructive or equational. 

The generic variables are provided as values of two attributes: 

parameters a list u of independent parameters 

vars a list x of dependent variables (equational proofs only) 

For equational proofs the variable lists x and u are chosen in such a way that u 
is a maximal independent set of variables. 

The basic attributes (with GeoCode values) are: 

Points the free points of the proof scheme 

coordinates assignments that compose step by step the generic geomet- 
ric configuration of the proof scheme 
conclusion the conclusion of the proof scheme 

This already completes the data required for a constructive proof scheme. For 
equational proof schemes the following additional attributes are defined: 

polynomials a list of GeoCode predicates that correspond to polynomial 
or rational conditions describing algebraic dependency re- 
lations in the given geometric configuration 
constraints a list of GeoCode predicates that correspond to polynomial 
non degeneracy conditions 

solution a way to solve the algebraic problem (given in extended 

GeoCode syntax) 

The proof idea can be sketched within the Proofldea attribute as plain text 
if not yet evident from the code. A CRef link to a PROBLEMS record with 
a UT^rjX formulation (kept in a separate place since one problem can be solved 
by several different proof schemes) and an optional Comment section are also 
available. 



3.3 The GeoCode Syntax 

The design of the generic GeoCode language is mainly inspired by our goal 
to fix proof schemes in a common format that easily translates to the special 
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input formats of different provers. It strictly distinguishes between identifiers 
(names for points, lines, circles) and symbols (names for the generic variables) 
to avoid the name space overlap that is typical for symbolic computations. We 
use Perl like syntax (i.e. , \$[a-zA-z] [a-zA-zO-9]* in Perl regexp notation) 
for identifiers and small letter / digit combinations (i.e., [a-z] [a-zO-9] * in Perl 
regexp notation - we don’t allow capital letters to avoid name clashes both in 
Reduce and Mathematica) for symbol names. 

At the moment we assume proof schemes to be composed by a sequence of 
assignments with nested function calls as right hand sides that refer to previ- 
ously defined geometric objects and scalars as arguments. This may change in 
the future in favor of a more XML-compliant syntax since several provers (and 
diagram drawing tools) don’t support nested function calls. 

Most CAS use parentheses both to group arithmetic expressions and argu- 
ment lists in function calls. Since this cannot be distinguished within a regular 
language we use the Mathematica convention (i.e., brackets) for function call no- 
tation 2 . The naming conventions support the identification of syntactical units 
via regular expression matching and their modification for the needs of different 
target systems (with special naming conventions, to avoid clashes with protected 
names as D,0,E etc.). Perl tools support such syntactical translations. 

The names and signatures of all the GeoCode functions are stored in the 
SymbolicData GeoCode table and can be extracted, extended and modified in 
the same way as other SymbolicData records. Two such GeoCode records are 
reproduced in figures 2 and 3. The first one corresponds to an ’inline’ function 
that requires a special implementation, the second one to a ’macro’ with a generic 
definition in GeoCode syntax as value of the code attribute that can be used 
to create an implementation automatically. For a complete description of all 
functions see the SymbolicData GeoCode documentation. 

Equational GEO records usually contain also a solution tag with a descrip- 
tion how the algebraic task can be solved. This description is fixed in an extended 
GeoCode syntax. Interface packages for Maple, MuPAD, Mathematica, and Re- 
duce to map these generic commands to the special CAS command syntax are 
part of the SymbolicData distribution. For details see [10]. 

3.4 The GeoProver Packages 

To run proof schemes written in GeoCode syntax with a geometry theorem 
prover requires an implementation of the GeoCode syntax for the target prover. 
We propose an approach in two steps: First, use Perl tools to translate the 
generic proof scheme into a syntactic form that is more appropriate for the target 
system (change square bracket, fix variable and function names, etc.). This is well 
supported by the SymbolicData actions concept and sample implementations 

2 Note that most of the arithmetic expressions were replaced by new geometric pred- 
icates eq_dist and eq_angle in the GEO records during preparation of version 1.3 
(finished after the ADG-02 conference) to emphasize the geometric nature of the 
proof schemes. 
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################################################################ 
# Record ’GeoCode/pp_line’ 



<Id> GeoCode/pp_line </Id> 

<Type> GeoCode </Type> 

<Key> pp_line </Key> 

<call> pp_line [$A :: Point , $B :: Point] :: Line </call> 

<verbose> line through A'and B </verbose> 

<description> 

The line through <math>A</math> and <math>B</math> . 
</description> 



# End of record ’GeoCode/pp_line’ 
################################################################ 



Fig. 2. The GeoCode “inline” record pp_line 



################################################################ 

# Record ’GeoCode/altitude’ 



<Id> GeoCode/altitude </Id> 

<Type> GeoCode </Type> 

<Key> altitude </Key> 

<call> altitude [$A :: Point , $B :: Point , $C :: Point] :: Line </call> 

<verbose> altitude from A'onto g(BC) </verbose> 

<code> ortho_line [$A,pp_line [$B, $C] ] </code> 

<description> 

The altitude from <math>A</math> onto <math>g(BC)</math> . 
</description> 



# End of record ’GeoCode/altitude’ 
################################################################ 

Fig. 3. The GeoCode “macro” record altitude 



for such translators are contained in the bin/GEQ directory of the SymbolicData 
distribution. 

Second, write an interface package in the language of the target prover that 
maps the GeoCode functions to the prover specific functions (and proving ap- 
proach). The author’s GeoProver packages implement such interfaces for proofs 
using the coordinate method and Maple, MuPAD, Mathematica, or Reduce. 

Below we give some examples of the interplay between GeoCode syntax, 
GEO proof scheme records, and their algebraic counterparts computed with 
the GeoProver package for one of the target CAS. For a formal description of 
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all functions we refer to the documentation [8]. For some target systems the 
GeoProver package provides also an undocumented plot extension that allows 
to draw diagrams from geometric configurations. 

A first prototype of the GeoProver grew out from a course of lectures for 
students of computer science on this topic held by the author at the Univ. of 
Leipzig in fall 1996. It was updated and completed to version 1.1 of a Reduce 
package after a similar lecture in spring 1998. Later on in cooperation with Malte 
Witte, at those times one of my students, the package was translated to the other 
target systems. 

Since version 1.2 there is a separate description of the GeoCode language that 
was fixed in SymbolicData format and added as GeoCode table to the Symbolic- 
Data project later on. Now the complete GeoProver source code is generated 
from a platform-specific ’inline’ code part for the basic functions and generic 
GeoCode code values for advanced functions using special SymbolicData tools. 
This facilitates a concise code management of the GeoProver implementations 
for different CAS if the GeoCode standard changes during development. 

4 Some Examples 

To demonstrate the efforts required to compile new GeoCode proof schemes, to 
translate them with SymbolicData tools to GeoProver applications and to run 
them on different target CAS we consider some examples. 

4.1 The “cathedral example”, [12, 5.3] 



P,Q are centers of the arcs CB,AC, 
respectively; arcs DE , DF are drawn 
with R, S as centers, respectively , and 
\AQ\ as the radius. Further, \AQ\ = 
\BP\ = ||AB|. The goal is to find the 
radius r of the circle tangent to arcs 
EA,ED,HM (with center A) and K M 
(with center D) as a function of |AB|. 

We take the problem formulation and notational conventions from [12] with 
s = \AB\ = 1. In that paper line AB is taken as z-axis with origin at M (written 
as $M in GeoCode syntax) and the following coordinates are assigned to points 
(using the GeoCode point constructor): 

$0 : =Point [0 , y3] ; $A : =Point [-3/ 12,0] ; 

$Q : =Point [7/12,0]; $M : =Point [0,0] ; 

Kapur’s original proof scheme contains 6 algebraic conditions that arise from the 
tangency conditions of two circle pairs. A circle tangency condition yields 3 poly- 
nomials: If Ti = (xi,y\) is the point of tangency of the circles around A and O 



c 




Fig. 4 : The cathedral example 
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these conditions are ’Ti on the circle around A\ ’Ti on the circle around O’, and 
’Ti, A, O are collinear’. 

Here is that statement in the formal GeoCode syntax: 

<vars> [xl,yl,x2,y2,y3,r] </vars> 

<Points> 

$0:=Point [0,y3] ; $A : =Point [-1/4 , 0] ; $Q : =Point [7/12 , 0] ; 

$M : =Point [0 , 0] ; $T1 : =Point [xl ,yl] ; $T2:=Point [x2,y2] ; 

</Points> 

<coordinates> 

$cl : =pc_circle [$A , $M] ; $c2 : =pc_circle [$Q , $A] ; 

</coordinates> 

<polynomials> 

$polys : =List [ 

on_circle [$T1 , $cl] , on_circle [$T2 , $c2] , 
is_collinear [$A, $T1 , $0] , is_collinear [$Q,$T2,$0] , 
r~2-sqrdist [$0 , $T1] , r~2-sqrdist [$0 , $T2] ] ; 

</polynomials> 

pc_circle is the center-point circle constructor. Note that on_circle [$T1 , $cl] 
and sqrdist [$T1 , $A] -sqrdist [$M, $A] yield the same algebraic translation. 
Hence a similar proof scheme may be composed without references to circle 
objects 3 . 

To solve the algebraic problem eliminate all variables but r from $polys 
and solve the remaining equation for r. This can be stored as the value of the 
solution attribute of the GEO record in extended GeoCode syntax in the fol- 
lowing way: 

<solution> 

$result : =geo_solve [ 

geo_eliminate [$polys , $vars ,List [xl ,yl ,x2 ,y2 ,y3] ] ,r] ; 
</solution> 

Less variables are required if Ti and Ti are taken as ’circle sliders’ that 
algebraically translate to rational parameterizations of such points: 

<vars> [xl,x2,y3,r] </vars> 

<Points> 

$0:=Point [0,y3] ; $A:=Point [-1/4,0] ; 

$Q : =Point [7/12,0]; $M : =Point [0,0] ; 

</Points> 

<coordinates> 

$T1 : =circle_slider [$A, $M,xl] ; $T2:=circle_slider [$Q,$A,x2] ; 
</coordinates> 

3 With GeoCode version 1.3 better use the (geometric) predicate eq_dist instead of 
the (algebraic) sqrdist difference. The former can be interpreted also, e.g., by a 
diagram drawing tool. 
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<polynomials> 

$polys : =List [ 

is_collinear [$A, $T1 ,$0] , r~2-sqrdist [$0 , $T1] , 
is_collinear [$Q , $T2 , $0] , r~2-sqrdist [$Q , $T2] ] ; 

</polynomials> 

Note that in this case the algebraic translations of the geometric conditions yields 
in fact rational functions since the coordinates of T\ and T% are not polynomial 
but rational. A special algebraic approach is required for such proof schemes 4 . 
A brute force call to solve ($polys , $vars) with MuPAD produces 24 solutions 
with 12 different values of r. 

A third approach uses the explicit circle tangency condition, that corresponds 
to a polynomial condition on the circle parameters. It requires only 3 variables 
and translates to a polynomial system. Take center O and a circumfere point X 
on the y-axis with generic y-coordinates 

$0 : =Point [0 , xl] ; $X : =Point [0 , x2] ; 

$A:=Point [-1/4,0] ; $Q : =Point [7/12 , 0] ; $M : =Point [0 , 0] ; 

add variables for the arcs (circles) that should be tangent 

$cl : =pc_circle [$A , $M] ; $c2 : =pc_circle [$Q , $A] ; 

$c3 : =pc_circle [$0 , $X] ; 

and fix the tangency conditions and another one for the radius r = x\ — X 2 of 
the circle C 3 

$polys : =List [ 

is_cc_tangent [$cl , $c3] , is_cc_tangent [$c2 , $c3] , r-(xl-x2)]; 

The problem is of equational type and poses a deduction task. There are no 
independent variables and an algebraic solution can be obtained if xi,X 2 are 
eliminated from the polynomials and the remaining equation is solved for r. The 
corresponding GEO record is given in figure 5. 

To translate and run that code with the GeoProver MuPAD version first call 
the SymbolicData MuPADCode action on the record Cathedral_l . sd 

symbolicdata MuPADCode Cathedral_l . sd 

It maps the proof scheme to MuPAD syntax, resolves name clashes and yields 
( GeoProver package loading omitted) 

//==> Example Cathedral_l 

clear_ndg() : 

delete ’xl 1 , ’x2’ , ’r’ ; 

_vars : =geoList (xl ,x2 ,r) ; 

// coordinates 

_0:=Point(0,xl) ; _X:=Point(0,x2) ; 

4 Rational expressions can be avoided using homogeneous point coordinates. 
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###################################################################### 

<Type> GEO </Type> 

<Key> Cathedral_l </Key> 

<prooftype> equational, deduction </prooftype> 

<vars> [xl,x2,r] </vars> 

<Points> 

$0 : =Point [0 , xl] ; $X : =Point [0 , x2] ; 

$A:=Point [-1/4,0] ; $Q : =Point [7/12, 0] ; $M: =Point [0 ,0] ; 

</Points> 

<coordinates> 

$cl : =pc_circle [$A, $M] ; $c2 : =pc_circle [$Q , $A] ; $c3 : =pc_circle [$0 , $X] ; 

</ coordinates> 

<polynomials> 

$polys : =List [is_cc_tangent [$cl , $c3] , is_cc_tangent [$c2, $c3] ,r- (xl-x2) ] ; 
</polynomials> 



<solution> 

$result : =geo_solve [geo_eliminate [$polys , $var s , List [xl , x2] ] , r] ; 

</ solution> 

###################################################################### 



Fig. 5. A proof scheme for the cathedral example [12, 5.3] as GEO record 



_A : =Point (-3/12 , 0) ; _Q : =Point (7/12 , 0) ; _M:=Point(0,0) ; 

_cl : =pc_circle (_A, _M) ; _c2 : =pc_circle (_Q , _A) ; 

_c3 : =pc_circle (_0 , _X) ; 

// polynomials 

_polys : =geoList (is_cc_tangent (_cl , _c3) , 

is_cc_tangent(_c2,_c3) , r-(xl-x2)) ; 

// solution 

_result : =geo_solve (geo_eliminate (_polys , _vars , geoList (xl ,x2) ) , r) ; 
quit ; 

Now we can run that script with MuPAD to get the solutions 
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similar to [12]. Note that r = — Al. corresponds to the position of X on the ’top’ 
of O since r = X\ — X2- Kapur missed the solutions r = They correspond 
to imaginary coordinates of O and hence are virtual. 

A similar computation yields the length of the radius of the circle in the top 
region of figure 4: With origin at D, the center 0\ and a circumfere point X of 
that circle on the y-axis we get the proof scheme 



<vars> 

<Points> 



[xl,x2,r] 



</vars> 
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$D : =Point [0 , 0] ; $01 :=Point [0,xl] ; $X:=Point [0,x2] ; 

$S : =Point [5/6 , 0] ; $P : =Point [-1/3 , 0] ; $B : =Point [1/2 , 0] ; 
</Points> 

<coordinates> 

$cl : =pc_circle [$S , $D] ; $c2 : =pc_circle [$P, $B] ; 

$c3 : =pc_circle [$01 , $X] ; 

</ coordinates> 

<polynomials> 

$polys : =List [is_cc_tangent [$cl , $c3] , 

is_cc_tangent [$c2,$c3] , r-(xl-x2)] ; 
</polynomials> 

Running the corresponding computation with MuPAD yields the result 
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4.2 The Generalized Steiner Theorem, [19, Ex. 7] 

Take three points A2, B2,C2 respectively on the three perpendicular bi- 
sectors of BC , AC, AB of any triangle ABC such that 



d{A 2 ,BC) = t ■ \BC\, d(B 2 , AC) = t ■ \AC\, d(C 2 ,AB) = t ■ \AB\, 



where d{P, QR) denotes the distance of the point P from the line QR 
and t is an arbitrary non-negative number. 

Then the three lines AA2, BB2, CC2 are concurrent. 




Note that the statement of a prob- 
lem may be included in a Symbolic- 
Data record as value of a new attribute, 
say Text, since each such record admits 
“undefined” attributes that are handled 
by the SymbolicData tools in the same 
way as “defined” ones. In the Symbolic- 
Data data base problem statements are 
stored in a special table PROBLEMS and 
cross referenced in GEO records since 
several proof schemes may refer to the 
same problem. In the same way biblio- 
graphical information can be related to 
GEO records. 



D.Wang solved that problem in [19] with oriented areas using a Clifford al- 
gebra approach that avoids the introduction of virtual solutions. A straightfor- 
ward coordinatization with independent variables u\,U2,t and dependent vari- 
ables X\,X2,X3 goes as follows: 
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<vars> [xl,x2,x3] </vars> 

<parameters> [ul, u2, t] </parameters> 

<Points> 

$A : =Point [0 , 0] ; $B : =Point [0 , 1] ; $C : =Point [ul ,u2] ; 

</Points> 

<coordinates> 

$A1 : =midpoint [$B , $C] ; $B1 : =midpoint [$A, $C] ; 

$C1 : =midpoint [$A , $B] ; 

$A2 : =line_slider [p_bisector [$B , $C] ,xl] ; 

$B2 : =line_slider [p_bisector [$A,$C] ,x2] ; 

$C2 : =line_slider [p_bisector [$A,$B] ,x3] ; 

</coordinates> 

<polynomials> 

$polys : =List [sqrdist [$A1 , $A2] -t~2*sqrdist [$B , $C] , 
sqrdist [$B1 , $B2] -t~2*sqrdist [$A, $C] , 
sqrdist [$C1 , $C2] -t~2*sqrdist [$A, $B] ] ; 

</polynomials> 

<conclusion> 

$con : =is_concurrent [ 

pp_line [$A,$A2] , pp_line [$B , $B2] , pp_line [$C, $C2] ] ; 
</conclusion> 

It yields a polynomial system with 8 solutions in {x\, X 2 , £3) in the rational 
function field k(u\,U 2 ,t) that correspond to the different combinations of ori- 
entations of the triangles ARC2, BCA 2 , CAB 2 - I checked this with Maple 8, 
Mathematica 4.1, MuPAD 2.5, and Reduce 3.7 and the solution 

<solution> 

$sol : =geo_solve [$polys , $vars] ; 

$result : =geo_simplif y [geo_eval [$con, $sol] ] ; 

</solution> 

geo_solve, geo_simplify and geo_eval are interfaces to the different syntax of 
the solve, simplify and subs commands of the target CAS. 

All four CAS found that for exactly two of the 8 solutions the theorem is 
valid (i.e., the $result simplifies to zero). Note that with t 2 replaced by t only 
Mathematica and Reduce found the same answer. Maple and MuPAD created 
expressions with several root symbols that could not be completely simplified 
in the following computation. Maple detected even reducible RootOf’s with the 
evala in its answer. 
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Fig. 7. The Miquel Circle on four con-cyclic points 



4.3 Miquel’s Axiom, [13, Ex. 5] and [4, Ex. 7.6] 

If four circles are arranged in a cyclic sequence, each two successive 
circles intersect at a pair of points, and a circle passes through one point 
in each pair, then the other four intersection points are also concyclic. 

Both [13] and [4] quote the problem as hard for the coordinate approach and 
propose a more geometric solution within their geometric framework. 

The problem has even a constructive solution: Take four points A ±, . . . , A 4 
on a circle around the origin O, centers Mi 2 , . . . . M 41 for circles passing 
through {Ai, A 2 ), . . . , (A 4 , A{) and compute for each consecutive pair of such 
circles the second intersection point (it has rational coordinates in the parame- 
ters Mi, U 2 , M 3 , Mi, ... , V 4 ). Here is the GeoCode formulation: 

<prooftype> constructive </prooftype> 

<parameters> [ul, u2 , u3, vl, v2, v3, v4] </parameters> 

<Points> $0 : =Point [0 , 0] ; $A1 : =Point [1 , 0] ; </Points> 
<coordinates> 

$A2 : =circle_slider [$0 , $A1 ,ul] ; 

$A3 : =circle_slider [$0 , $A1 ,u2] ; 

$A4 : =circle_slider [$0 , $A1 ,u3] ; 

$M12 : =line_slider [p_bisector [$A1 , $A2] ,vl] ; 

$M23 : =line_slider [p_bisector [$A2 , $A3] ,v2] ; 

$M34 : =line_slider [p_bisector [$A3 , $A4] ,v3] ; 

$M41 : =line_slider [p_bisector [$A4 , $A1] ,v4] ; 

$cl2 : =pc_circle [$M12 ,$A1] ; 

$c23 : =pc_circle [$M23 , $A2] ; 

$c34:=pc_circle [$M34,$A3] ; 

$c41 : =pc_circle [$M41 ,$A4] ; 

$B1 : =other_cc_point [$A1 ,$cl2,$c41] ; 
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$B2 : =other_cc_point [$A2 ,$cl2 , $c23] ; 

$B3 : =other_cc_point [$A3 , $c23 , $c34] ; 

$B4 : =other_cc_point [$A4,$c34, $c41] ; 

</coordinates> 

<conclusion> 

$result : =is_concyclic [$B1 , $B2, $B3, $B4] ; 

</conclusion> 

The simplification of the resulting rational expression in seven parameters indeed 
turned out to be very hard and only Maple 8 (about 4 sec.) and Mathema- 
tica 4.1 (about 250 sec.) mastered the task 5 . The situation completely changes 
if Vi, . . . , V 4 are assigned random integers (not too big, < 100). Several runs with 
different settings always yield 0 also with MuPAD or Reduce. 

An even harder challenge is the Miquel-5 example [4, Ex. 7.6] that starts with 
five points forming a pentagon. There is a straightforward constructive proof 
scheme but only Maple 8 (about 40 sec.) could simplify the resulting rational 
expression in reasonable time. 



5 Implementing the Geo Code Standard 

For real usability of the Geo Code concept one has to estimate the efforts re- 
quired to implement this standard for a given prover. Even though there is not 
yet practical experience with already existing geometry theorem provers the se- 
mantic similarity of “foreign” special geometric proof schemes suggests that such 
an interface could easily be implemented. We suggest to divide that implemen- 
tation in two parts as described above: The first part provides (e.g., Perl based) 
tools to translate GEO proof schemes into a prover specific form that fixes re- 
quirements of naming and syntax conventions. In a second part these translated 
proof schemes are passed to a special interface that maps the Geo Code to the 
prover’s functionality. 

In general, for the second part one has to implement the GeoCode function- 
ality in the target prover language. This requires extensibility of that language 
and access to the source code or interaction with the system developers if such an 
interface cannot be added as a supplementary package. Let’s consider the proof 
scheme language described in [3] as target and estimate the efforts required to 
construct such a supplementary interface package to the GeoCode language on 
top of it. 

First, note that the GeoCode syntax provides not only points but also line 
and circle objects as geometric primitives. It is a special design decision of [3] 
(and other geometry theorem provers, e.g., D. Wang’s GEOTHER [16]) not to 
introduce the latter objects as basic but to take only point objects as primi- 
tives. The interface package implements lines and circles as pairs of points 6 with 

5 On a 800 MHz Pentium PC with 256 MB RAM under Linux. 

6 To support such an approach and also to make the GeoCode language more “geo- 
metric” we removed direct constructors for lines and circles from their homogeneous 



